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A COMPUTER SOLUTION FOR THE DYNAMIC LOAD, 
LUBRICANT FILM THICKNESS AND SURFACE 
TEMPERATURES IN SPIRAL BEVEL GEARS 


SUMMARY 

A complete analysis of spiral bevel gear sets is presented. 
The gear profile is described by the movements of the cutting tools. 
The contact patterns of the rigid body gears are investigated. The 
tooth dynamic force is studied by combining the effects of variable 
teeth meshing stiffness, speed, damping and bearing stiffness. The 
lubrication performance is also accomplished by including the effects 
of the lubricant viscosity, ambient temperature and gear speed. A 
set of numerical results is also presented. 
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CHAPTER I 
INTRODUCTION 

Spiral bevel gears, found in many machine tools, automobile 
differential gears, and aircraft gears, are important elements for power 
transmissions. However, the basic mechanisms which govern the major 
failure modes of spiral bevel gears are not fully understood. Unlike 
spur or helical gears, the complicated geometry of spiral bevel gears 
makes the problem considerably more difficult. In military applications 
such as gear transmissions used in recent V/STOL aircraft, gears 
are often designed under conditions very close to the failure limits 
to gain a maximum power density. A thorough understanding of spiral 
bevel gears under critical operation is urgently needed to prevent 
premature failure* 

Basically, there are two classifications of gear 
failure, one of which is structural failure, which includes flexure 
fatigue, tooth breakage and case crushing, and the other .is the 
lubrication related failure which includes wear, surface pitting and 
scuffing. 

Structural failures are usually related to the material properties 
and improper geometrical design or unexpected overloading. While the relia- 
bility of spiral bevel gears can be improved by better material 
selection and better geometrical design, it can also be enhanced by 
a more accurate understanding of the dynamic loading between a pair of 
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teeth as it travels through the contact. 

Since the basic lubrication process in bevel gears is not clearly 
understood, the lubrication related failures .particularly in surface 
pitting and scuffing, are much more difficult to predict. Current methods 
used in industries to predict lubrication related failures are not 
satisfactory because they are mainly empirical. It has been shown in 
failure tests of gears and rollers that the surface pitting and 
scuffing are affected critically by lubricant film thickness and surface 
temperature. A tool that can predict film thickness and surface tempera- 
ture in the gear teeth contact is welcome in design against pitting and 
scuffing. 

This paper presents an analysis to predict the dynamic loads 
between gear teeth contact which has a contact ratio greater than unity 
and a variable stiffness along the path of contact. In addition, an 
analysis is presented to predict the minimum film thickness and 
the surface temperature under given dynamic loads. It is based on the 
most recent theories on film thickness and traction in elastohydrodynamic 


contacts. 
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CHAPTER II 

TOOTH PROFILE AND CONTACT PATTERN 
OF SPIRAL BEVEL GEARS 

2.1 Introduction 

Recent developments in the finite element method (FEM) and 
elastohydrodynamic lubrication (EHD) have stimulated a growing interest 
in the basic research of the dynamics and lubrication process in spur 
and helical gears (1-6). This interest is now gradually being extended 
to spiral bevel gears for which there seems to be little known about 
the tooth load sharing and lubrication characteristics, especially at 
high speeds and heavy loads. 

Unlike involute spur and helical gears whose surface geometry 
is relatively simple and readily derivable, the surface geometry of 
spiral bevel gears is extremely complex owing to its three dimensional 
characteristics. For studying the performance and failure of spiral 
bevel gears, it is essential first to find an efficient method to 
describe the surface geometry of a pair of mating gears in terms of the 
tooth form and the tooth contact patterns. Although this method is 
currently available at Gleason Gear Works (7-11), it is not fully suffi- 
cient for studying the dynamic load and lubrication analysis • It is the 
intent of this research to extend the existing surface geometry analysis, 
known as the Tooth Contact Analysis (TCA) , to generate additional 
quantities needed for the lubrication and dynamic analysis. 
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The formulae for obtaining the coordinates of tooth surfaces 
based on the Gleason generating system a re derived in Section 2.2. The 
position of each point on the surface is described in cylindrical 
coordinates in terms of its position along the rotational axis of the 
pinion or gear, its radial distance from the axis, and its angular posi- 
tion about the axis from an arbitrarily chosen reference plane. The 
equations are solved by an iterative method. 

The contact points of a pair of mating gears and pinions for a 
running position are determined in Section 2.3. The sliding velocity 
and the direction of the unit normal vector at each contact point can 
be obtained. In addition, the principal radii of the mating surfaces which 
characterize the geometry at the contact point can be computed. 

Some numerical results of a typical case are shown in Section 
2.4. The effect of different running positions on contact patterns 
are studied. Finite element grids of a typical spiral bevel tooth 
profile are presented. 
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Determination of Tooth Surfaces for 

Generated Spiral Bevel Gears 

In Tooth Contact Analysis (TCA) the gear and pinion tooth 
surfaces are mathematically described by the machine settings and the 
cutter specifications used in making the gears. Before proceeding to 
describe TCA, it is first desirable to review some basic concepts of 
bevel gearing. The principal reference planes of a spiral bevel gearing 
system are shown in Fig. 2-1. The axial plane of a single gear is the 
plane containing the gear axis and a given point on the tooth surface. 
The pitch plane is the plane perpendicular to the axial plane and tan- 
gent to the pitch circle. The transverse plane is perpendicular to the 
axial plane and the pitch plane. Figure 2-2 and Figure 2r-3 show, 
respectively, the sectional views in the axial and transverse plane of 

a bevel and pinion system, and illustrate the most common terms used in 
bevel gears. 

In order to understand the formulae derived in this chapter, a 
brief explanation of the cutting process is in order. Figure 2-4 
illustrates the Gleason cutting process — 1) the wheel cutter used in 
cutting spiral gears, 2) the cradle which positions the cutter in 
reference to the work and carries it through its generating motion, and 
3) the blank into which the teeth are cut. The dresser, which is used 
to restore the cutting action of the wheel by fracturing and tearing 
away the dull grains to expose fresh cutting edges or clear away the 
imbedded material, is shown in Figure 2-5. Before the teeth are cut, 
the wheel cutter is predressed to shape the cutter. The relation 
between rotation of the cradle and rotation of the work is controlled 
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by the rotation of a cam which rotates about a fixed center. 

The formulae derived here are based on the kinematics of two 
contacting bodies. There is no relative velocity in the normal direction 
at the cutting point and the teeth are assumed to be rigid. The 
formulae are written in vector notation. Most of these vectors are 
self-explanatory except for the following two terms: 

I(g,e) R , a rotation of vector A through an angle 0 about the 

unit vector g. 

(-,"t,n) T , a 3 x 3 matrix to transform a vector into a new 

coordinate system having the base vectors indicated 
in the parentheses. Since the right hand cartesian 
coordinate is used , it is only necessary to indicate 
two base vectors. 

2.2.1 Basic Formulae 

To describe the manufacturing process, one needs to know the 
relationship among the rotating elements, and the relative movement 
between the work and these elements. 

2.2.2 Data Defining the Gear Tooth Surface 

The gear tooth surfaces are controlled by the machine settings. 

The symbols of data required are listed in Table I. Sample data are 
given in Appendix A. These values are usually tabulated on a Gleason 
Gear Engineering Standard form. 

A spiral bevel gear and pinion are always of opposite hand. The 
formulae considered here apply to a set consisting of a left-hand pinion showi 
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TABLE I 


SYMBOLS OF DATA FOR MACHINE SETTING 



Gear 

Pinion 

Concave 

Convex 

Number of teeth 

N 

n 


Face width 


F 


Shaft angle 


z 


Outer cone distance 


A. 




0 


Addendum 

a 




s 

d 

P 


Dedendum 

b 

b 



S 

P 


Pitch angle 

r 

T 


Face angle 

r 

T 



0 

0 


Root angle 

r 

x 



r 

r 


Spiral angle 


*M 


Nominal wheel diameter 


D 




0) 


Machine center to back 

X 

X 

X 


g 

P 

P 

Sliding base 

V 

X bp 


Blank offset 

E 

E 



mg 

mp 

I j 

mp 

Cam setting 

T 

T 

T 


g 

P 

P 

Eccentric angle 

e 

s 

6 p 


g 

P P 

Cradle angle 

Q g 

% 

% 

Cam guide angle 

a 

a 

ct 


g 

P 

P 

Standard cam setting 

T n 

T 

T 


go 

PO 

PO 

Cam pitch radius 

Ou 

g 

nu 

P 

nu 

p 

Index interval 

n . 

n. 

n 


i g 

IP 

ip 
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“ Fl8ure 2-6 and right-hand gear shorn in Figure 2-7. The drive side 
is assumed to be. the concave side of the pinion running with the convex 
side of the gear. The opposite combination is the coast side. Since 
the gear tooth data define the side under consideration, these general 
formulae apply to either the drive side or the coast side. 

2.2.3 Mechanism of Tam 

The relative rotation of the work and cradle is controlled by 
the index head and cam. The index head is used to transfer accurately 
the rotation ratio between the cam and workpiece. By describing the 
motion of the cam with a specific index head setting, one can obtain 
the rotations of the work and cradle, and their angular speed ratio. 

The follower of the heart-shaped cam is fixed in the cradle, and 
the cam rotation center is fixed. Figure 2-8 shows the position for 
a standard cam setting. As the cam rotates about its own rotation 
center, the fixed follower forces the cradle to rotate about its center. 
The required rotation of the work is obtained by using the proper 
indexing. The relative motion can be readily changed by using a 
different cam geometry, cam setting, cam guide angle and index interval. 
The symbols used in the motion formulae are: 
e cam rotation 

* work rotation 

V cradle rotation 

° cam guide angle 

An u change in cam setting 


cam pitch radius 
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T cam setting 

X standard cam setting 

o 

For a standard heart-shaped cam, the mechanism of motion of the 
cam and cradle at center of roll can be thought of as a pair of internal 
gears which has a pitch radius n u for the pinion and DC + n u for the 
gear (see Figure 2-9). If the cam rotates an angle e Q , the cradle will 

rotate an angle n u /(n u + DC) times e q « 

The change in cam setting An u is defined as the actual cam 
setting T minus the standard cam setting T q (see Figure 2-10) . The 
dashed lines are for the standard cam setting for a rotation e Q . For a 
change in An u> the following relations are derived. 


+o = DC* 


u 


( 2 . 1 ) 


n e 
u o 


An 

sin if) =■ sin sin(e Q - # Q ) 


( 2 . 2 ) 


where i|r and g are the actual rotations for the current setting. 

A center line of the follower is parallel to and is m the middle 
of the follower surfaces. The cam guide angle is defined as the angle 
between this center line and the line from the cradle center to the 
cam rotation center before the cam rotates (see Figure 2-lla) . Because 

the cradle, the distance from the cradle center 


the follower is fixed in 
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to the center line of the follower i 


s constant and is always equal to DC 


sin a. 


Figure 2-llb shows the relative positions with a caw guide angle. 
Using e o as the standard cam rotation angle, one obtains 


ib . = — e 

v 0 DC +n e 0 
u 


(2.4) 


sin 1^2 - sin. ip^ + sin a 


(2.5) 


f " a 


( 2 . 6 ) 


e 9 = - tfi 

2 0 v 2 V Q 


(2.7) 


e = e 2 “ a 


( 2 . 8 ) 


If both the cam setting and the cam guide angle are changed, one 
can find the following relations based on the geometry shown in Fig. 2-12 


for a standard cam rotation s : 

o 


ib - u e 

0 DC + n 0 

u 

An 

(2.9) 

sin \l> 2 = sin \p Q + sin 

“ + "5c sin(e o - V 

(2.10) 

^ - a 

* z 


(2.11) 

2 o y 2 *2 


(2.12) 

e = e _ - rt 

2 


(2.13) 
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X 



(2.14) 


The ratio between the cradle and the work rotation becomes 


= _n_ 5* (2.15) 

d\ 


where 

N M = n u COS *0 + An u cos(e 0 " V 


(2.16) 


°M = N M + DC COS ^2 

These ratios of the rotational velocity and the rotation of the work and the 
cradle are necessary for computing the velocity and the normal direction 
at the contact point. Their derivations are shown in Appendix B. 

2.2.4 Basic Machine Set-up and Vector 
Coordinate System 

The cradle settings in the spiral bevel gear cutting machines 
of the Gleason Works are shown in Fig. 2-13. The settings are fixed during the 
cutting process. From this figure, the following relations are obtained. 


s = 16 sin 
P 2 


q po 



Q + 360° 
P 


6 

s = 16 sin 
g 2 


q g0 



(2.18) 

(2.19) 

( 2 . 20 ) 


( 2 . 21 ) 
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Vectors in these formulae have been taken with respect to a 
set of right-hand rectangular cartesian coordinates, x, y, z fixed in 

the machine plane viewed from the front of the cutter toward the gear 
as shown in Fig. 2-14. 

In order to be consistent with the vector formulae, the follow- 
ing matrix for unit vectors, d, e, and c is used. 

d = (1, 0, 0) 

e = (0, 1, 0) (2.22) 

c = (0, 0, 1) 

2.2.5 Point on the Gear 

In cutting the gear, the moving parts include the rotating work, 
cradle, wheel cutter, cam, and dresser. The work rotation, cradle 
rotation and cam rotation are related to each other. It is desirable 
to choose cam rotation as the independent variable. Thus, there are 
only three independent rotations which define a specific point on the 
gear tooth surface, i.e. the cam, dresser, and wheel cutter. 

The procedure to determine a point on the gear tooth surface is 
to choose arbitrary values for the cam rotation and wheel cutter 
rotation, then find a value of dresser rotation such that the kinematic 
requirement of contact is satisfied, i.e., the relative velocity at the 
cutting point is perpendicular to the normal of the point being cut in 
the gear tooth. 

Symbols shown in Fig. 2—14' for all position vectors are 
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defined as: 

A : vector from the machine center to the point on the surface 

5 • vector from the machine center to the intersection of plane 

containing the cutter tips and the wheel cutter axis. 

R : vector from the crossing point to the point on the surface. 

B : vector from wheel center to the point on the surface. 

7 : vector from the intersection of plane containing the cutter tips 

and the wheel cutter axis to cutter tips. 

n : unit normal vector to the cutter blade. 

t : unit vector along the direction of dressing diamond movement, 

g : unit vector along the gear axis. 

"p : unit vector along the pinion axis. 

The initial position of these vectors when q g = 0, 0 g = 0 are 
g = (-cos r r , 0 , -sin r.) 


Y = “Y j e 
'go dg 


S = S d 
go g 


n * ( 0 , cos <(» , sin <J> ) 
go g g 


t = ( 0 , sin <{> , -cos 4 > ) 
go g g 

Because the convex side of the gear tooth is formed by an inside 
blade of the cutter, by using the sign convention given below, the same 
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formulae (used in computer programs) can be used for both sides: 


Drive Side 


+g ' £ *d g ' % 


Coast Side 
♦ - - ♦ 


dg 


X, = X, . 
dg dig 


X, = - X. 
dg dog 


L , = L, . 

dg dig 


L dg L dog 


Reasonable initial values of cam rotation e„ (standard cam 

Og 

setting) and wheel rotation 0 are 

g 


e o g ■ 0 


( 2 . 23 ) 


0 - q 0 — (1; 

g go 


( 2 . 24 ) 


where iji^ is the spiral angle, using e from cam formulae to obtain 

^g’^g and ratio of cradle rotation to work rotation. The new position 
vectors become 


q g = q go + \ 


S = S -(c, q V 
g gO v ’ 4 g' 


r =* r (c, q -0 ) 
g go' H g g ; 


‘is ■ v c ’ vv 


R 
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n = n(c, q -6 ) R 
lg g° ’ g 


d i g ’ V c - VV 


B „ = r - (0.2 - a . ) c + X. n- - (2.0 - E ) d 
gO g dg' dg lg dg lg 


A g0 S g + B g0 


R g0 = X g g - Emg e - X feg c + A gQ 


where B n : vector from wheel center to dresser center. 

g° 

A n : vector from machine center to dresser center. 

g° 

R A : vector from crossing point to dresser center. 

go 


After setting certain values for the two independent variables 

e and 0 , one can determine the dresser angle 0 , that satisfies the 
0g g 

kinematic requirement of contact between tooth surface and cutter. 

To find this dresser angle 0^ , an initial value (0^ ~ 0) fi rst 
assumed. From Fig. 2-14, the new position vectors are 


d 2g “ d lg (n ig’ 9 dg ) 


C 2G “ ‘ig^lg' V 


A = A + L, d 
g gO dg 2g 


R = R n + L , d 0 
g gO dg 2g 
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B = B n + L , d. 
g gO dg 2g 


The normal of the cutter blade at this instant (given * , 0 values) is 

Og eg 


n = {cxB)xt } 
g g 2g unit 


(2.27) 


The velocity of the cutt er relative to the work at the point being cut 


is 


c x A g + g x R g (2.28) 

g 

The product n g • V is checked to see whether or not it is zero; if not, 

the secant method is used to iterate 0 dg until n g • V = 0 is satisfied 

for the prescribed values of e rt and 0 . 

°g g 

The tooth surface can be described completely with the formulae 
derived so far by varying the values of e Qg and 6 g independently. 

2.2.6 Point on the Pinion 


V = it 

g cA j 


The procedure to find the pinion tooth surface is exactly the same as 
that for the gear tooth surface, except for some initial vectors which 
are different. From Fig. 2-15, the mathematical forms are described 
briefly as follows : 


Drive Side 

Coast Side 

♦p ' V 

* P ■ v - *+ 

% ■ V 

X - - X.. 
dp dip 

^dp ^dop 

dp dip 
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The initial position vectors are: 


p = (-cos y r , 0, -sin yj 


S A - S d 
p0 P 


r = r , e 
p0 dp 


n = (0, -cos <)> , -sin <j> ) 
pO P P 


t . = (0, -sin 4, cos 4 ) 
pO P P 


Assume e ft - 0, 6 = q - i|» ; one can calculate X ,* and ratio of 

Op p pU M P P 

cradle rotation to work rotation by cam formulae described earlier. 


S ■ So - S 


S = s(c, -q ) 

p pO p 


r ■ r (c, 0 -q ) 
p pO P P 


t. - t (c, 0 -q ) 
lp pO P P 


n = n (c, 0 -q ) 
lp pO P P 


d i P * d(c> S ■\> 


B = 7 - (0.2 - z. ) c - X. n - (2.0 - E ) d 
po p dp dp lp dp lp 


A n - s + B 
pO p p0 


R n - x p + E e - X. c + A 

p0 p D P 
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change 0 , until V • n =0. 
a P p 

The tooth surface of pinion is thus determined by changing the 
values of and 0^ independently. 
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2.3 Tooth Contact Behavior and Motion Graph 

The only tooth curvature of involute gears with an infinite 
lengthwise radius of curvature is the tooth profile curvature, so the 
instantaneous tooth contact pattern will be a straight line extending 
from one end of the tooth to the other. If elastic deformation is 
included, the contact pattern will be an elongated rectangle, from one 
end of the tooth to the other. If the teeth are crowned, lengthwise 
curvature is introduced. This results in a point contact for rigid 
teeth, or an elliptical contact for elastically deformed teeth. 

Theoretical conjugate bevel gear teeth can be manufactured, but 
they are sensitive to shaft misalignment and housing deformation. Due to 
this, bevel gears are produced with a lengthwise and a profilewise mis- 
match. Therefore, the instantaneous tooth contact pattern will be a 
theoretical point for rigid teeth or an ellipse when elastic deformation 
is included. 

By choosing the proper amount of mismatch in both directions for 
a particular case, edge contact and excessive contact pressure can be 
avoided in the presence of misalignment due to shaft and housing 
deformations. 

The purpose of the motion graph is to describe the nature of the 
motion transmitted and the effects of contact between adjacent pairs of 

teeth on the tooth contact pattern. 

The gear contact problem is essentially a Hertzian contact prob- 
lem, with the pressure distribution in the Hertzian ellipse a main concern 
in EHD. The shape of the Hertzian ellipse can be readily determined 
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from the principal radii which can be computed analytically from the 
gear teeth geometry. 

2.3.1 Matching a Contact Point 

A method to find the contact point by machine setting and r unn ing 
position is derived in this section. 

In addition to the rectangular cartesian coordinate used -in 
Chapter II, a point on the tooth surface with respect to the tooth axis 
and the crossing point can be specified in the axial plane by: 

L : the distance along the tooth axis from the crossing point. 

R i the redius from the tooth axis . 

The primary effect of misalignment is to shift the tooth contact 
to the boundaries of the tooth surface, resulting in edge contact and 
extreme loading conditions. Figure 2-16 shows a pair of 
spiral bevel gears; the following errors in relative running position of 
the pinion with respect to gear are discussed: 

2 r ' running shaft angle. 

Xpr : pinion apex withdraw -which is the linear displacement in 
the H-direction (along the pinion axis) . 

^g r : 8 ear a P ex withdraw -which is the linear displacement in the 
^“direction (along the gear axis) . 

E r : offset-which is the linear displacement in the V-direction. 
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Consider a point on the gear tooth surface obtained by assumed 

values of e and 9 as described in Section 2.2, in order to match 
Og g 

the gear axis from its manufacturing position to a proper running 
position relative to the pinion axis • Suppose one lets the manufacturing 
axis of the pinion coincide with the pinion axis at running position. 

Then the gear axis at running position is as shown in Figure 2-17 , and 
the new position vectors relative to the new gear axis are: 



(2.32) 


(2.33) 


where y is the desired rotation angle of the gear from a specified 
S 

reference. In the computer program, n is set at zero to find the first 

o 

contact point, and then is increased to obtain the successive contact 
points. 


Assuming that a point P on the gear surface obtained by a pair 
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of values and 0^ is a contact point under running conditions, the 
position vector for this point expressed in the pinion axis can be 

obtained by placing the gear in running position relative to the pinion 
as shown in Figure 2-18 



= X p 
P r 


+ E e - X g. +R 
r gr 1 gl 


(2.34) 


One can define this contact point P in the pinion without con- 
sidering its angular position by 




(2.35) 


In the manufacturing procedure, there is a set of particular 

values of e Qp and 0 p to produce a point in the pinion with the same 

values of L and R as those of L and R . 

P P pi pi 


The mathematical way to find these particular values of e and 

Op 


0 p is described as follows. 

Assuming arbitrarily e Qp and 9^ and using the method described 
in section 2.2, one can get R p which is the vector from crossing point 
to the contact point Q between cutter blade and work. The vector R 


can be specified by 


P 


L - 
P 


R = 
P 



* P 
x p | 


(2.36) 
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Iterate and 9 simultaneously until L - L and R - R ., that is, 
Op p P P-L P P-L 

the same distance around the pinion axis from crossing point is obtained 

for point Q and point P. 

Then rotate the pinion axis to make point Q coincide with the 
contact point P; the rotation angle required, An pJ is shown in Figure 
2 - 19 . 


A rip = sin 



(R . x p) 
P 1 


( 2 . 37 ) 


R 


Point P is obtained under the assumption that arbitrary values 

of e and 9 were made. Under a given tooth profile and running posi- 
Og g 

tion, not every point in the tooth could be a common contact point. A 
common contact point must have the same normal direction for both pinion 
and gear. 

Assuming that a point P is the common contact point, the normal 
direction at this point on the gear for and 9 g is on the 

pinion it is 

n = n (p,An ) R (2 ' 3 

Pi P P 


where n is obtained at point Q. n , and n are compared to see if 
p g-L P-L 

they agree to within a specified error. If not, the secant method is 

used to to change £ and 9 simultaneously to find another contact point 

og g 

until the two normal directions agree. This procedure yields a point of 
contact between gear and pinion. 
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In fully matched gears, the constant velocity ratio is main- 
tained. In single mismatch case (lengthwise mismatch) the theoretical 
contact point is a point of the original surface, so the uniform velocity 
will be maintained. There is only one theoretical contact point in 
combined mismatch, the point of the original surface; therefore, the 
velocity ratio will not remain constant, and there is a theoretical 
instantaneous velocity change at each changeover. 

The motion graph shows the error in angular displacement of the driven 
member against rotation of the driver. It is assumed throughout that 

the pinion is the driving member and that it is rotating at a constant 
velocity . 

2.3.2 Condition of Contact 

The surface topography of spiral bevel gears cannot be readily 
expressed mathematically by an analytical function. However, one can 
use second-order surface approximations to describe the surface around a 
point of contact. By this approximation the principal radii and relative 
surface velocity at the contact point can be analytically determined. 
Approximating a surface to the second order gives 

2 

z = 0.5 ax + b x y + 0.5 c y^ (2.39) 

The numeric coefficients were chosen' so that 


a 


curvature in x-direction 
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b = 


3x9y 


twisting 


(2.40) 


c 



curvature in y-direction 


express this surface equation in compact form (s) 



Except for the tooth breakage, the tooth failure modes are 
mainly due to local contact pressure which causes surface pitting or 
scoring at high sliding velocities. Therefore, to develop reliable 
formulae for gear design it is necessary to have knowledge of contact 
pressure and surface velocity on the gear tooth. 


2. 3. 2.1 Surface Coefficients 

Figure 1-20 shows the surface at the point of contact for the gear. 
The tangent plane coordinate system is chosen in such a way that the 
plane containing the normal direction vector of contact point (n g ) an< ^ 
the vector c of wheel axis is the yz plane of the tangent plane coordinate 
system as shown. Let t coincide with y axis and n g with z axis, so 

o 

the vector T x n is in the direction of the x-axis. 

§ s 

In this chosen tangent plane coordinate system, there is no 
twisting for the surface of the generator which is the wheel surface 
around the point of contact (b gwg = 0) and the curvature in the x-axis 
direction of this generator surface (a gug ) can be derived based on the 
geometry shown in Figure 2-20. 


sin ip = n _ * c 
mg g 


(2.41) 
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r „„ = I B x c | sec 1 1> 
ng g 1 r oig 


scog r 


ng 


2 - 21 . 


The curvature in the y-axis direction (c ) is shown i 

swg 


t - n tan i jt — c sec ib 

g g wg <Og 


sin i - n • (t x n ) 
g lg 2g g 


sin o - t, • (t x n ) 
“g 2 g g g 


x = L sin a 0 + 

dg cog 


y = L cos a 0 + 

dg wg 


A = 0.5 L (sin i )0 2 - 0.5x 2 — 
a & g r 

ng 


2 L 

— 0.5 9 [-L. sin i - — ^ sin^ a ] 
dg g r 


ng 


0)g J 


but 


4 - 0.5 y 2 ji- 

0)g 


hence 


1-1 Sin \ 1.2 


P L. 
wg dg cos a 


-r -2 — tan' a 

2. r ojg 

cog “8 


(2.42) 

(2.43) 
Fig. 

(2.44) 

(2.45) 

(2.46) 

(2.47) 


(2.48) 
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(2.49) 


swg p 


where o is the radius of curvature in y-axis direction, 

wg 

In the actual running position, the z-axis of the tangent plane 
is n , and the y-axis is 

gl 


t „ = t (g,p) R (e, AZ) R 
gt g 


(2.50) 


For the pinion, the surface of the generator around the point 
of contact is shown in Fig. 2-22, and surface coefficients a g , b SU)p 
and c can then be obtained: 

SU)p 


sin * up - -n p • c 


(2.51) 


t = n tan t|> „ + c sec ip 

p p dip Wp 


(2.52) 


sin i = n. • (t_ x n ) 
p lp 2p p 


(2.53) 


sin a = t„ * (t x n ) 
ajp 2p P P 


(2.54) 


r = |B x c sec ip _ 
np p 1 ^P 


(2.55) 


, , sm i 1 , 

_i l - — 2- - JL tan 2 a 

p L, 2 r up 

H u>p dp cos a np 


(2.56) 
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a 

scop 


1 


r 

np 


b = 0 

sup 


( 2 . 57 ) 


= 

^ P 0) P 


In the real running position, the tangent plane coordinate 
in the z-axis direction is 


n 


= n 


Pi gl 


( 2 . 58 ) 


and in the y-axis direction is 


V ' V p ’ 4 V 


( 2 . 59 ) 


One uses the vectorial velocity and acceleration of a known 
point on the wheel and angular velocity vector of the wheel about that 
point to describe the movement of the generator. This point is chosen as 
a point of contact at a given instant. The motion parameters are: 


V : linear velocity vector, 

a : linear acceleration vector. 


u 


angular velocity vector. 
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The above are relative motion vectors denoting the motion of a moving 

generator relative to a generatee. 

The values of these motion parameters can be computed for the 
gear as follows (see Fig. 2-23 for direction of vectors): 




d 2 t] 

0) = 
eg 

W* ’ 

O 

dX 2 


(2.60) 


oo = „„ c + 8 

g eg 


(2.61) 


E = R - A 
g g g 


(2.62) 


V = g x E + « g x A g 


(2.63) 


3g = g x (g x E g ) + (o» cg g x c + a c§ c) x A g 


( 2 . 64 ) 


+ a) x (cu x A ) 
g 8 8 


(2.64) 


Transforming these vectors to tangent plane coordinates 4 


_ — — X 

CD = u) (- , t , n ) 
tg g ’ g 8 


v tg ■ V-- V v 1 


(2.65) 


a tG + a G (_ * V n g )T 
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Similarly, for the pinion, one obtains 



0) = 0) C + D 

P cp 


E = R - A 
P P P 


V = p x E +a) xA 
P P P P 


a p * P x (p x E p ) + (u)^ p x c + a cp c) x Ap 


+ m x (u» x A ) 
P P P 


( 2 . 66 ) 


(2.67) 


( 2 . 68 ) 


(2.69) 


(2.70) 


Transforming to the tangent plane coordinate system, 


“ tn = “ (“. t , n ) T 
tP P P P 


V. = V (-, t , n ) T 
tp p P p 


(2.71) 


a - a (-, t , n ) T 
tp p p p 


The problem is, then, to determine the three components of the 
^^ erence surface matrix (Aa, Ab, Ac) in terms of the given surface 
components and the motion parameters. From 2.12, the following relations are 
obtained for the gear (subscripts x,y,z denote the component of a vector): 
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tan a 


id. + a V 
tgy sa)g tgx 


se iii - c V 

tgx S0)g tgy 


( 2 . 72 ) 


a + V id - V id 
tgz tgx tgy tgy tgx 


u - c V 
tgx s(Dg tgy 


( 2 . 73 ) 


(D - C V 

Ar „ tgy so)g tgy 

g v + V. - V tan a 

6 g tgy tgx sg 


( 2 . 74 ) 


Ab = - Ac tan a 
g g sg 


( 2 . 75 ) 


Aa = - Ab tan a 

g g sg 


( 2 . 76 ) 


The coefficients of the surface of the gear around the contact point are: 


a = a + Aa 
sg s«Jg g 


b = Ab 
sg g 


( 2 . 77 ) 


c = c + Ac 
sg suig g 


for the pinion: 


tan a 


id + a V 
tpy soap tpx 


sp (D - C V 
v tpx s<Dp tpy 


( 2 . 78 ) 
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v 

P 


a. + V u — V m 
= tpx tpy tpy ^tpx 

0) - c v 

tpx stop tpy 


( 2 . 79 ) 


AC P 


v 

P 


to 


tpx 


+ V 

tpy 


- c v 

stop tpy 

- V tan ot 
tpx sp 


Ab = -Ac tan a 
P p sp 


Aa = -Ab tan a 
P P sp 


( 2 . 80 ) 


The coefficients of the surface of the pinion in the pinion tangent plane 
coordinate system are 


a 


sp 


a 

scop 


4- Aa 

P 


b = Ab 
sp p 


( 2 . 81 ) 


c 

sp 


c + Ac 

Sli)p p 


It is necessary to rotate the pinion tangent plane coordinate system to 
coincide with the gear tangent plane coordinate system in order to express 
the coefficients of both surfaces in the same coordinate system. The angle 
required for rotation is A which can be found from 


sin A = t 

Pt 





( 2 . 82 ) 


and then to transform the coefficients of the pinion surface into this 
coordinate system, 
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2 2 

a ■ a cos A + 2 b cos A sin A + c sin A 

spg sp sp sp 

2 2 

b = (c - a ) cos A sin A + b (cos A - sin A) (2.83) 

spg sp sp' sp 

2 2 

c = c cos A - 2 b cos A sin A + sin A 

spg sp sp sp 

2. 3. 2. 2 Geometric Contact Ellipse 

Using a value of 0.00025" for the separation at the extremities 

of the contact and the surface curvatures of each member at the point of 
contact, the size and direction of the contact ellipse are determined 
as (Fig. 2-24). 



Ab = b - b 

spg sg 


(2.84) 


Ac 


c — c 
spg Sg 


tan 2g 


2Ab 

Aa - Ac 


The two values of g are the principal directions. Assuming that g 2 *- s 
first or fourth quadrant, the curvatures in the g^ and g 2 directions are: 
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— « Aa cos g 1 + 2 Ab cos ^ sin + Ac sin 2 


(2.85) 


12 o 

— = Aa cos g 2 + 2 Ab cos g 2 sin g £ + Ac sin 


P]_ anc * P 2 are P r i nc ip a l radii of the difference surface of gear and 
pinion surface, the lengths of the semi-axes of the contact ellipse in the 
tangent plane are: 



( 2 . 86 ) 


2. 3. 2. 3 Velocity of Contact Point 

Consider that the gear surface is stationary. During a small unit 
time, the contact point on pinion moves from P to P' with velocity V 
due to rotation of the pinion. During this same time interval, the 
point of contact moves on the pinion surface to P" with a velocity v 
(Fig. 2-25), the velocity of contact point is the sum of these two 
parts. 


u = V + v 
s 


(2.87) 
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The pinion point velocity relative to the gear for unit pinion 
rotational velocity is 


V = m.(R x g) + (R x p) 
si 1 gl pi 


Transforming this vector into the common tangent plane 


v s - V-' v v T 


The angular velocity of the pinion body relative to the gear is 


« 1 - -P - g 


Transforming this vector into the tangent plane, 


“ * "i ( -> V v )T 


The surface coefficients of gear and pinion in this tangent plane are 


(s ) - 
g 


a b 
sg sg 


b c 
sg sg 


( a b 

spg spg 

b c 

spg spg 


Since p" is to be a point of contact, the rate of change of 
the normals of the gear surface and the pinion must be the same: 
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_ / V 

dt u ' ,S g' f° r the gear 

= - v(s p ) + a) x n for the pinion 

where u x n is due to the rotation of the pinion surface, and -v(s ) is due 

P 

to curvature of pinion surface. Hence, 

-u(s ) = - (u - V.) (s ) + « x n 

© op 


u[(s ) - (s)]=uxn + V (s) 
8 P s p 


u(As) = co x n + V (s) 
s P 


u = [oixn + V (s)] (As) 
s p 


( 2 . 93 ) 
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2.4 Results and Discussion 

One typical set of spiral gears currently being tested at NASA 
Lewis Research Center was used to illustrate the numerical calculations. 

The machine settings, which are provided by Gleason Works, are shown in 
Appendix A. The points on the tooth surface were obtained by procedures 
described in Section 2.2. The contact condition was determined by 
the method shown in Section 2.3. A series of solutions were obtained 
to simulate the spiral bevel gears subjected to a range of operating 
conditions, and results are presented here. In order to obtain a 
unified view, all the curves are plotted in the axial plane. 

2.4.1 Tooth Profile 

The methods in Section 2.2, which use actual cutting motions 
of the gear cutting tools, result in an exact description of the tooth- 
forms as functions of spiral angle, blade-edge radius, etc. 

A typical tooth profile for the gear shown in Fig. 2-26 was 
plotted in Fig. 2-27. For the purpose of applying FEM in gear design, 
a typical finite element model for a single gear tooth was shown in 
Fig. 2-28. 

2.4.2 Contact Pattern 

The spiral bevel gears are relieved both profilewise and length- 
wise, so that there is only a single contact point in the center of the tooth 
This combined mismatch permits the gear pair to tolerate displacements 
under load and assembly errors. Theoretical point contact is, of course, 
only an imaginary concept. Under even light loads, tooth surfaces 
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compress and deflect to cause surface contact. This surface contact 
area is roughly elliptical. The sum of those elliptical contact areas 
along the path of contact is called "Tooth Bearing." Correctly made 
teeth will show a tooth bearing centrally located on the tooth surface, 
when the gears are in their standard relative running positions. The 
amount of lengthwise mismatch is readily observed by the length of the 
tooth bearing. Profile mismatch is not observed as readily owing to 
the fact that the tooth height is relatively small; nevertheless, the 
existence of profile mismatch is evidenced by the lack of any heavy 
concentration of bearing at the tips of the teeth. 

The results illustrated here were obtained by analyzing the 
effects of running position errors on the path of contact. The errors 
discussed here include pinion apex withdraw (X^) , gear apex withdraw 
( X gr ) » runni ng shaft angle (^ r ) and running offset (E^) . Their effects 
are calculated by changing the value of a single term while keeping the 
others unchanged, and these results are shown from Fig. 2-29 to Fig. 

Th e path of contact moves to the toe as the pinion apex withdraw 

increases and moves to the heel as it decreases. The shape of the path of 

contact remains the same, but is inclined more as it approaches both 

sides. The gear apex withdraw will produce the same effect, but in the 

opposite direction. As increases, the path of contact moves to the 

heel, and the opposite occurs as X decreases. Both of these 

8 r 

have as their chief effect a pressure angle error; the tooth contact 
will move toward the tooth tips on the gear tooth and toward the tooth 
flanks on the pinion. When the shaft angle becomes larger, the contact on 
both the gear and the pinion shifts to the toe. At the same time, the contact 
moves toward the tooth tips. The contact goes the opposite way when the 
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shaft angle decreases. This effect is expected. The principal effect 
of the offset is to tilt the tooth in its mating slot, so that the con- 
tact shifts to the toe on one side and to the heel on the other side. 

Usually, the running position errors are introduced simul- 
taneously; the results on the drive side are investigated for the sliding 
velocity and radii of the difference surface at the contact point with 
three separate sets of errors which result in the center contact, heel 
contact, and toe contact. The tooth bearing for these three contact 
cases is plotted respectively on the same graph. Assuming the combined 
Hertzian deformation to be 0.00025 M , which corresponds approximately to 
to the case of a pair of lightly loaded gears, and^using the curvatures 
of the contacting surfaces, the size and direction of the contact 
ellipse are determined. The larger the radius, the bigger the ellipse. 
The sliding velocity has also been plotted at selected points. See 
Fig. 2-34 and Fig. 2-35. 

2.4.3 Motion Graphs 

The rotation errors of the driven member against the driving 
member for the entire duration of contact for several adjacent meshing 
pairs were plotted. The corresponding velocity errors were shown in 
the same graph. Fig. 2-36 shows the curves on the drive side with a 
running position which results in a center contact. Fig. 2-37 shows 
the results on the coast side with the same running position. Figs. 

2-38 and 2-39 illustrate the contact patterns for the motion curves 
shown in Figs. 2-36 and 2-37. 
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The velocity curves in those figures do not intersect; at the 
transfer points there is a sudden jump in the velocity with which the 
driving tooth strikes its mate. This value is, therefore, a measure 
of the impact occurring at each tooth mesh. As mismatch is increased, 
the impact increases in magnitude, resulting in noise and stress. It 
is therefore essential that mismatch be kept to the smallest value 
compatible with the adjustability to the running position errors. This 
concept of motion errors is only applied to the case under a very light 

load. The elastic deformation due to moderate load will smooth out these 
errors • 
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pig, 2-1 . Principal Reference Planes 




Fig. 2-2. Sectional View in The Axial Plane For A Pair 
of Spiral Bevel Gears 



Fis. 2-3. Sectional View in The Transverse Plane For A Pair 
Of Spiral Bevel Gears. 
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Fig- 2-5. Movement of the Dresser 
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Fig. 2-7. Spiral Bevel Gears 




47 



Fig. 2-8. Standard Cara Setting 



Fig. 2-9. The Movement of Standard Cam Setting- 



Fig. 2-10. The Movement of Cam With a Change in Cam Setting 



( a ) 



C b ) 


/am With a Cam Guide Angle oi 
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Cradle Center 



p^g < 2-12". The Movement of Cam With Change in Setting 
and Cam Guide Angle <* 


5 
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3. Cradle Settings in 
Machine of the Glea 
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Fig. 2-21 . Wheel Surface Curvature (Pinion) in y Direction 
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Fig. 2-25. Surface Velocity 
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2-21, Tooth Profile Plotted by Calculated Points 
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Fig. 2-28. The FEM Model For a Single Tooth 
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Rota cion 

Fig. 2-36. Motion Graphs (drive side) For Successive Tooth Meshes 







75 


CHAPTER III 

DYNAMICS OF SPIRAL BEVEL GEARS 

3.1 Introduction 

One of the main uncertainties in gear failure analysis is the 
load imposed on the gear teeth at high speeds where the inertia forces 
of the gear wheel and shaft become significant. The dynamic response of 
the operating system plays an important role in gear design because the 
movement of the gear body will shift the contact bearing to an undesirable 
position as a result of the misalignment of shaft. The dynamic load in 
gearing has received continuous attention in the past, and many studies 
have been reported (1-9). As far as the dynamic response of bevel gears 
is concerned, most of the studies have been conducted experimentally in the 
laboratory. One of the recent measurements is by Terauchi and Fujii (10) 
in which the dynamic load of gear teeth of straight bevel gears is 
obtained by using a power circulating gear testing machine. There is 
a lack of analytical study of the vibration of spiral bevel gears. For 
instance, some industrial practice still depends upon Buckingham's formulae 
or upon empirical factors for designing spiral bevel gears. 

Analytical studies of spiral- bevel gears are extremely complex 
because: 1. the tooth profile of spiral bevel gears is generated by the cutting 

machine, and it is difficult to describe the tooth surface, which is much 
more complicated than that of spur and helical gears; 2. the tooth contact 
patterns of spiral bevel gears vary with the assembly errors caused by 
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the apex error, shaft alignment, etc.; 3. being different from spur 
gears, the spiral bevel gears have also lateral vibration in addition 
to rotational vibration. In this report, efforts are made to explore 
the effects of all these factors on the dynamic behavior of the spiral 
bevel gears. 
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3.2 Formulation 

The main objective in this section is to formulate a dyn am ic 
model and equations of motion to simulate, in the steady state, the periodic 
motion of both the pinion and the gear as well as the tooth load during 
a typical cycle, in which a pair of teeth traverses through the zone of 
action from point A to point C as shown in Fig. 3—1. 

In developing the equations of motion, the pinion and the gear were 
assumed to be rigid bodies each having six degrees of freedom. The 
supporting bearing was assumed flexible with known spring stiffness. At 
the contact of each pair of teeth in the zone of action, the teeth were 
assumed to be connected by a linear spring which is oriented normal to 
the contact point and has a stiffness to be determined separately by a 
finite element model. 

The equations of motion of a rigid body in three dimension are 
(see Fig. 3-2): 


where 


Z F = m a 


(3.1) 


£ *<G * H G 

F : external force vector 
m : mass 

a : acceleration vector of the mass center 
Mq : vector of external moment about the gravity center 


Hq : vector of rate of change of angular momentum about 
the mass center 
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The following derivation is described separately for pinion and gear. 


3.2.1 Pinion 

Let the moving coordinate axes x p , y p , z^ be the principal axes 

of the inertia of the body so that the axis of the rotation of the body 
will direct along one of these axes, which is named the dominant axis. 
The pinion shaft is not fixed in the moving coordinate, but rotates 
about the dominant axis y p at an average velocity u p which is the 

operational rotation velocity of the shaft. XYZ axes are fixed co- 
ordinates in space (see Fig. 3—1) . In the case of the pinion, the 
dominant axis is along the Y axis, hence the angular momentum of the 
pinion about the mass center is 



I , 9 , i' + I , (6 , + u> )V 

x p x'p p y p y p p p 


+ I , e , k' 

z p z'p p 


(3.2) 


and the rate of change, of H_ in the moving coordinate 


H = 1 e i' + i, e j' + I , 0 k' 
Gp x p x'p p y P y p p z P z’p p 


(3.3) 


+ a x (i 0 , i' + i , (0 , + u )j' + i » 0 . k ') 

p x'p x'p P y p y'p p p z P z P p 


8 


(0 , i' 
x'p p 


+ 0 


+ e_. k') 


, 1 ' T V I 

y P P z p p 


(3.4) 



Cl x'p 0 x'p 


— I I U> 0 t ■*" 

y p p z'p 


(I , - I , )0 , 0 » J 
v z'p y'p' y P 2 P 


i' 

P 
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•[i , 0 , 

[_ y p y 


+ (I , - I , )0 , 9 , 

p x p z'p x p z'p 


jp < 3 - 5 > 


•fi , e , 

[_ z P z ? 


„ + i , <d e , + (i , -i , )e , e 

p y p p x'p y ' p x'p x'p 


• -i 

9 , 

y p 


k f 

P 


®yip is the rotation velocity variation from to . The true 

* • 

rotation velocity along the dominant axis is 0 , + to . Assuming 0 , 

y P P y p 

• • • 

is much less than co^ and the magnitude of 0 x » p > 9 y i p and 0 , is of the 

same order, and neglecting the higher order terms, one obtains 
a linearized equation 


H 


Gp 


I , 0 i — X , to 0 ■ 
x'p x'p y'p p z'p 


i' 

P 


1,9, 

yp y p 


(3.6) 


1,6, + I . to 0 , 

z p z p y'p p x'p 


k' 

P 


so the governing equations become 


M N 

I * -r + l F . = m (x' I' + y' j' + z’ k') 
i=l rpl i=i C P 1 P C P P cp J p cp p / 


(3.7) 


M _ _ N _ _ 

1 Y_ . xF .+ y Y .xF + Input Torque 

i=l rpx rpi ^ cpx cpi 


-[ 


1,0, - I , to 0 , 

X P x p y'p p z'p 


] 


i' 

P 


(3.8) 
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+ [*1,0,1 j' 

L y p y p J p 

(3.8) 

+|i , e , + i , uj e , 1 k' 

[_ z’p z'p y'p p x'pj p 


where H is number of bearing reaction force, N is number of contact 

force, and x' , y* , and z' are the components of the absolute 
cp cp cp 

acceleration of the pinion mass center in the x' , y' and z' direction. 


3.2.2 Gear 

These results are similar to the equations for the pinion, except 
that the dominant axis in the case of the gear is X axis and the operational 
rotation velocity is id , the equations of motion for the gear are: 


M 


l 

i=l 


F . + 
rgi 


N 

l 

i=l 


cgi 


= m 


(x* i* + y' j' + z' k’) 

r* rr " /-> /> r> o rr 


g eg g 


cg-'g 


eg g 


(3.9) 


M 

l 


i=l 




+ Output Torque 


i’ 

g 


+ ^ I y ' g% ' g + I x’ g u g 9 z' g] * 8 

+ J^ I z , g 9 z , g I x , g ai g 9 y , g] 8 


(3.10) 
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Some preliminary work must be completed before the twelve 
equations in the system can be solved. These are; 1) expressing the 
force terms F^ and in terms of the displacement unknowns in the 
fixed coordinate; 2) expressing the displacement unknowns in the 
fixed coordinate by the displacement unknowns in x’y'z* system in which 
the equations are solved. 

Because the contact force is assumed to always be in the normal 
direction of the contact point on the tooth surface, the scalar quantity of the 
contact force can be used effectively in the derivation if the direction 
vector of the contact point is known. For simplification, the matrix 
notation is used, and some basic column matrices are introduced: 


{D'} = (x' , y', z', 0 , , 6 , 0 , ) T 

g 8 g g x'g’ y g z'g' 


{D'} = (x', y', z' , 0 , , 0 , , 0 , )■ 

P p p’ x'p’ y p z'p' 


(3.11) 

(3.12) 


which are the displacement unknowns in the moving coordinate system, 


while 


{ D } = (x ,y ,z ,0 ,0 ,0 ) 

g g g g xg’ yg’ zg' 


(3.13) 


{ D } « (x ,y ,z ,0 9 ,0 )' 

P P P P xp, yp zp 


(3.14) 


are the displacement components in the fixed coordinate. 

From the configuration of the gear system, the displacement 

components at the contact point i in the normal direction (G ) can be 

n 
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expressed in terms of the rigid body displacements as: 


{G .} = [DG.] {D} 

ni J g i J g g 


(3.15) 


(G ,} - [DG ] {D} 

ni p i p P 


(3.16) 


and the bearing reaction force at bearing j can also be expressed 
in terms of rigid body displacements by: 


{F . } = -[DK. ] {D} 

r 3 g j S g 


(3.17) 


{F . } = - [DK.] {D} (3.18) 

r 3 p IP P 

where i and j are from one to the number of contact forces and number of 
bearings, respectively. The displacements at contact point i due to 
elastic deformation along the direction of contact normal are: 


{6 .} = [DC 1 {F .} 

nci g ij g cj g 


(3.19) 


(6 . } = [DC,.] {F .} 

nci p ij p cj p 


(3.20) 


where i, j are from one to the total number of contact points. The 
unit outward normal vector at contact point i, is denoted nci. The 
positive value of F ^ is taken so that contact force is acting along 
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the same direction as the outward pointing normal vector. Since only the 
scalar quantity of the contact force is considered in the equation, the 
relation {F ci > = {F ci } p is always true at each contact point. 

The total displacements in the normal direction, which include 
both the rigid body and the elastic deformations, will be the same for 
the two contacting bodies. This relation gives: 


n . * + {G.} ) - - n . 

cl g i g i g exp 


({6.} + {G. } ) 

x p x p 


(3.21) 


+ {G n-T } „ = . } — (G } 

nci g ni g nci p ni p 


[DC..] (F . } + {G . } - - [DC .] {F } - {G } 

Jig ci g ni g 1 jiV cip ni p 


{F } = - ( [DC . . ] + [DC..] ) 1 ({G .} + {G . } ) 

C1 g Ji g Ji P nx g nx p 


(3.22) 


(3.23) 


and 


{F ci>g - - ‘[’Vg + WjiV 


-1 


([DG.] {D} + [DG.] {D} ) (3.24) 

1 g g i P P 


which is the equation relating the contact force and the rigid body 
displacement unknowns. The displacements in different coordinate 
systems can be transferred from one to another by means of the following 
relations (10) , 
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( 3 . 25 ) 


( 3 . 26 ) 


( 3 . 27 ) 


where [0 U ] is the transformation matrix for rotation in x’y' 2 ' coordinate, 

k 

[0 ] is the transformation matrix for rotation in xyz coordinate. 

Incorporating all the above force and displacement relations 
into the equations of motion for the pinion and gear and then rearranging 
the matrices, one can obtain the final form of the twelve equations. 


[m] 





{R} 


( 3 . 28 ) 


The details of matrix [m] , [c] and [k] are derived in Appendix C. 

After examining the equations of motion, it was found that the 
fourth equation governing the rotational motion along the dominant axis 


85 


of the gear, and the eleventh equation governing the rotational motion 
along the dominant axis of the pinion can be reduced to a single equation . 

The reduced equation is similar to that describing the vibratory system of 
two masses connected by a spring. In this equation, it is convenient to convert 
the rotational motion into a motion along an artificial line of action. 

Like the line of action in spur gears, the artificial line of action in 
spiral bevel gears is defined as the effective component of y x ii to 
produce the torque along the dominant axis, where y is a position vector 
from the mass center to the contact point, and n is a unit normal 
vector at the contact point* 

The 4th and 11th equations of motion are rewritten here as 

VgVg = Y cgi X ^x'g + (Output Torque) x , g (3.29) 


I , 0 , 

y p y p 


K Y C pi x F cpiVp + (In P ut T ° r 1“> y .p 


(3.30) 


The following new variables are introduced: 


0 = RBG • 0 , 

8 x'g 


(3.31) 


0 = RBP 

P 


y'p 


(3.32) 


Letting FMA denote the magnitude of the contact force and T be the static 

c 

force, the equations become 
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(I Y . x F .) , = FMA • RBG 

cgi cgi x g 


(Z Y . x F .) . = - FMA • RBP 

'cpi cpi y p 


(3.34) 


(Output Torque) , - 


(3.35) 


(Input Torque)^, p - T c • RBP 


(3.36) 


RBG and RBP are treated as the radii of the effective base circle of spiral 
bevel gears. Eqs. (3.29) and (3.30) become 


I x . tt •• 

9 * RBG • FMA - T • RBG 

RBG g c 


(3.37) 


I , .. 

-Zl£. e = 

RBP p 


- RBP • FMA + T * RBP 
c 


(3.38) 


Defining 


J = Ix ' g 
Xg RBG 2 


(3.39) 


J - 
yP RBP 2 


(3.40) 


one obtains 


J 0 = FMA - T 

xg g c 


(3.41) 


J 0 = - FMA + T 


(3.42) 
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Introducing an equivalent mass J and a relative displacement 9 

Pg 


J • J 

j = _*g ZE 

J + J 

xg yp 


(3.43) 


9=0-9 

Pg P g (3.44) 


the reduced equation, including damping terms, becomes 

J 0 + C 0 + FMA = T 

Pg pg c 


(3.45) 


After the relationship between the 11 displacement unknowns and the 
magnitude of the contact force is rewritten, the system can be reduced to 
11 second order differential equations with state dependent coefficients. 
These equations were solved numerically by the Runge-Kutta method. The 
unknown initial conditions (11 displacements and 11 velocities) can be 
iterated by taking the previous calculated values at the end of each 
iteration. The criterion of the iteration is to force the equation to 
satisfy the periodic contact condition. The static displacements and 
zero velocities are taken as first guessed values. 

To implement the above equations, the following state coeffi- 
cients are needed: 

(1) The tooth contact position as a function of the relative 
rigid body displacements of the two shafts, 

(2) The direction of the normal vector at the contact point, 

(3) The combined stiffness of the teeth 


at the contact point. 
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The teeth contact position and the direction of the normal vector at the 
contact point can be obtained from Chapter II. Because of the geometric 

complexity of spiral bevel gears, there is no satisfactory 
formula available to calculate the deformation of the tooth surface. 

A finite element program, to be described in Section 3.3, was used in 
this study to calculate the deformation due to a unit load applied at a 
given contact point of spiral bevel gears. 
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3* 3 Tooth Deflection 

For most gear applications, the contact ratio is greater than 
one, that is, there will be more than one pair of teeth in contact 
during some portion of the whole engagement. When the load is shared 
by two contacts, it cannot be assumed that the load is distributed 
equally among the pairs of teeth in contact because this is a statically 
indeterminate case. Therefore, one must consider the tooth deflection 
under the load for each pair in order to determine the load sharing 
characteristics among the pairs. 

Because of the complexity of the spiral bevel gear geometry, 
there are no existing simplified methods for calculating the suitable 
tooth deflection. In order to investigate the system response, the gear 
shaft must also be included in the calculation of tooth deflection. 

Some of the recent applications of finite elements in determining 
tooth deflection can be found in Refs. 9 and 12, where it was shown that 
more accurate results can be obtained by using the finite element method. 
Most of this work dealt with two dimensional problems and did not 
include the whole gear body. Fig. 3-4 shows a typical 8-node solid 
element grid pattern for a gear and a pinion with three adjacent teeth 
attached to the gear wheel and shaft. Fig. 3-5 shows the central 
tooth and its attached ring element of the gear. Fig. 3-6 shows parts 
of the gear shaft and the gear wheei. Fig. 3-7 shows the entire ring 
element with three adjacent teeth of the pinion. Fig. 3-8 shows the 
elements of the pinion shaft. The central tooth was subject to load to 
calculate the deflection. There are 941 nodes, 562 elements for the 
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gear model, and 1029 nodes, 584 elements for the pinion model. Using 
this grid as the input, one can readily compute the deflection 6 under 
a load P applied at any grid point in the tooth surface. For this 
analysis, the MARC-CDC program was used, the boundaries were considered 
to be fixed for all the points connected to thrust bearing to eliminate 
rigid body displacement, and the boundaries nodes connected to the 
radial bearing was allowed to move in the direction of dominant axis. 

The stiffness at grid point i is defined as 



(3.46) 




(3.47) 


The stiffness (KS , KS p ) of the point other than the grid point in the 
tooth surface can be calculated by the interpolation method. The 
details are shown in Appendix D. The combined stiffness at contact 
point is expressed as 


KS 


KS 


KS, 


• KS 

E 

+ KS 

P 


(3.48) 
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3 * 4 Results and Discussion 

A series of solutions were obtained to simulate the dynamic 
response of a set of spiral bevel gears currently being tested at NASA- 
Lewis. The data for this gear set and the lubricant data are listed 
in Table 3.1. Effects studied include running speed, shaft misalignment, 
and system damping. These results are presented in this section. The 
dynamic response is expressed by a dynamic load factor defined as the 
ratio of the maximum dynamic load along the contact path to the average 
static load. This factor is plotted as a function of speed with differ- 
ent damping ratios and contact ratios. 

3.4.1 Dynamic Load Variation 

For a constant input torque, the load on the contact point of 
the two meshing teeth along the path of contact is not constant. This 
load variation is mainly caused by the following factors: 

1. The variation of stiffness along the contact path. 

2. The transition from a single pair of contacts to a double 
and from double to single. 

3. The effective radius is not constant along the contact path. 
Fig. 3-9 shows the variation of effective radius of pinion. Fig. 3-10 

shows the variation and jump of stiffness for the transition of 
contacts. 

The main excitation to the gear system comes from the periodical 
change m teeth stiffness due to the alternating engagement of single 
and double pairs of teeth. The frequency of this excitation force, 
expressed as a meshing frequency, depends on the operating speed. 
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Therefore, it dominates the resulting mode of vibration. Fig. 3-11 to 
Fig. 3-14 show dynamic load variation at four different speeds in the 
case of central contact, that is, when the contact path is located 
centrally between the toe and the heel of the tooth. 

Since there are eleven degrees of freedom in the system, eleven 
resonating frequencies of the system should exist. In the low speed 
region where the excitation frequency from the change of stiffness is 
much lower than all resonating frequencies, the dynamic load response 
along the path of contact is somewhat like static load superimposed by 
an oscillatory load due to the resonating frequency of the system. 

When the speed is near the resonance region (Fig. 3.11), the 
dynamic load response becomes very severe (Fig. 3-12 and Fig. 3-13). 

The maximum dynamic load is much higher than the static load, which is 
the case when overloading occurs. Sometimes the oscillation of dynamic 
load will make meshing teeth separate when the load becomes negative, 

and thus generates noise and surface fatigue. 

As the speed increases beyond the zone of resonating frequencies, 
the dynamic load becomes smoother along the contact course, and the 
value is less than the static load (if the contact ratio is above one). 
The variation of dynamic load in this region is out of phase with the 
change of the teeth stiffness (Fig 3-14) . 

3.4.2 The Effect of Shaft Misalignment 

When the assembly errors are introduced in the system, the 

either end of the tooth surface (11) . 


contact bearing shifted to 
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^8* 3-15 shows the typical paths of central contact, toe contact, and 
heel contact. Usually the central contact is desired because it can 
tolerate more possible running position errors and avoid edge contact. 
The dynamic load response of toe and heel contactsare shown in Fig. 

3-16 and Fig. 3-17. The change of the contact bearing from the center 
to either edge will also change the contact ratio of the system because 
the tooth surface is not a perfect involute along the profilewise 
direction and is mismatched along the lengthwise direction. In the 
current example, the contact ratio for the toe contact is 1.26, 
central contact 1.16, and heel contact 1.0. In this case, if the con- 
tact bearing is moved farther toward the heel region, there would be no 
tooth contact between the time when the previous tooth finishes the 
contact and current tooth goes into the contact zone (discontinuity 
in tooth mesh) . This situation would cause a very large impact force 
which would generate noise and severe damage to the tooth surface. 

The effect of the tooth contact ratio on dynamic response is shown later. 

3.4.3 Contact Path Variation Due to Dynamic Response 

In addition to showing the contact paths due to assembly errors 
in the system in Fig. 3-15, the real contact path, including the effects 
of both assembly errors and running position errors induced by the 
dynamic responses, is plotted in the same figure. When 
this real contact path is compared with that caused by the assembly 
errors and running position errors induced by the average static elastic 
deformations, the deviation is found to be surprisingly small. One 
explanation of this small difference might be that the displacements 
change due to the dynamic oscillation are small, and they do not produce 
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a large change in contact path compared to those caused by the static 
displacement only. The closeness between these two contact paths suggests 
that one can use the average static elastic deformation to calculate the 
contact path, which can then be used directly to solve for the dynamic load 
and lubrication behavior without having to solve the dynamic load 
and contact path simultaneously using an iterative technique. The 
elimination of this iterative procedure greatly reduces the computation 
time. 

3.4.4 Effect of Speed 

Once the physical conditions of a gear set are determined, the 
dynamic response depends on the operating speed. For a system with one 
degree of freedom, such as spur gears, the maximum dynamic load occurs 
when the meshing frequency, which depends on the operating speed, is 
near the system natural frequency. Some peaks in the dynamic load 
are caused by the varied meshing stiffness along the contact path, and 
they appear at meshing frequencies lower than the system natural fre- 
quency. The dynamic load factor, defined as the ratio of maximum 
dynamic load to the average static load, is plotted against the gear 
speed to illustrate the effect of speed in Fig. 3-18. Since there are 
eleven degrees of freedom in the spiral bevel gear system, more peaks 
in the dynamic load are expected. 

It is seen that the highest dynamic load appears to occur near 
the natural frequencies which correspond to the mode associated with 
the largest displacements in the motion along the line of action. The 
frequencies marked t in Fig. 3-18 shows the system natural frequency 


c 


-a- 
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causing a larger displacement in the motion along the line of action and 
the ones marked 4- shows the system natural frequency with a small 
displacement in that motion. It is clearly shown that the dynamic load 
factor at the frequencies marked + has a peak response, and the response 
at the natural frequencies marked + is not necessarily a peak. 

3.4.5 Effect of Contact Ratio 

The contact ratio is defined as the ratio of the time required for 
one tooth to go through the whole contact zone to the time required for a 
periodic meshing cycle. It is believed that the load sharing 
characteristics caused by more than one tooth in contact will reduce the 
static load. The variations in the dynamic load factor due to the effect 
of changing contact ratio is shown in Figs. 3-19(a) to (c). It can be 
seen that the maximum dynamic load factor does not change much. However, 
the contact ratio*s effect is significant in high speed regions, where 
the load is spread out between meshing teeth pairs. A typical dynamic 
load variation with a high contact ratio along the contact path is shown 
in Fig. 3-20. 


3.4.6 Effect of Damping 

Since the damping forces are usually not known in the gear system, 
three arbitrary values are chosen for the damping coefficients: 2627, 

4378, and 6129 N sec/m (15, 25, and 35 lb.sec/in). These values are 
selected to give a range of non-dimensional damping ratios corresponding 
to those commonly used in spur gears (0.1-0. 2). The non-dimensional 
damping ratios corresponding to the above damping coefficients are 0.087, 
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0.14 and 0.203. The dynamic response for these damping cases can be 
observed from Figs. 3-18, 3-19(a), and 3-21. It is expected that the 
larger the damping force is , the smaller the dynamic load factor will be 
in the resonance region. The large damping force will also level off the 
peak of dynamic load factor in the subresonance region, and there is no 
effect on the dynamic load factor due to the damping force in the 
superresonance region. 

3.4.7 Effect of Bearing Stiffness 

It is well known that the bearing stiffness plays an important 
part in the dynamic load response because it directly controls the static 
displacements which determine the contact path. A large stiffness for 
supporting bearings is sometimes desirable because it pushes all the 
resonant frequencies beyond the range of the operating speeds. In Fig. 
3-22, the dynamic load factor is calculated for a system with "infinite” 
bearing stiffness. The contact ratio for this case is 1.16, and the 
damping ratios are 0.058 and 0.14. Since the bearing stiffnesses are 
infinite, only the rotational mode prevails. Three peaks are shown to 
exist in Fig. 3-22; two of these are the subharmonics. These results for 
infinite stiffness are compared with the results shown in Fig. 3— 19(a) 
for a finite stiffness of 3.5 x 10® N/m (2 x 10^ lb/in.). It is seen 
that the natural frequency of the rotational mode with an infinite 
stiffness is increased to 290 rad. /sec. from the natural frequency of 200 

g 

rad. /sec. for a finite stiffness of 3.5 x 10 N/m. 
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CONTACT PATH 


Fig. 3-1. Dynamic Modeling and Zone of Action of Bevel Gears. 
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Fig. 3-4. A Typical Section of Elements of Gear 
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Central Tooth Elements And Attached Kira Elements of 
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Pig. 3-7. The Elements of Rim And Three Adjacent Teeth of 
Pinion 
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5-10. Stiffness Variation Along The Contact path 
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Fife. 3 - 13 . Dynamic Load Variation Along The Contact Path at 
Gear Speed 300 rad/aec For Central Contact. 
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Fig. 3-14. Dynamic Load Variation Along The Contact Path at 
Gear Speed 523 rad/aec For Central Contact. 
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Dash Line For Contact Path Without Elastic Defor- 
mation;Solid Line for Contact Path With Elastic 
Deformation 
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Speed 523 rad/sec For Toe Contact 
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Fig. 3 - 17 . Dynamic Load Variation Along The Contact Path at 
Gear Speed 523 rad/sec For Heel Contact 
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Fig. 3-18. Dynamic Load Factor is Plotted Against Gear Speed 
With All Damping Coefficients = 15 lbf • s/in 
Contact Ratio = 1.16 
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Fig. 3-19. (b) Dynamic load Factor is Plotted Against Gear 
Speed With Contact Ratio =1.33 
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Fig. 3-21. Dynamic Load Factor is Plotted Against Gear 

Speed With All Damping Coefficients » 35 lbf • s/in 
and Contact Ratio = 1.16 
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CHAPTER IV 


LUBRICATION OF SPIRAL BEVEL GEARS 


4.1 Introduction 

One of the most important failure modes in lubricated machine 
elements is scuffing, which is where local wear occurs between the 
contacting elements when the lubricant fails to separate them and the 
protective thin surface film breaks down due to the high temperature. In 
gearing systems this failure is affected mainly by gear geometry, speed, 
load, and lubrication. For many years, the methods used to predict the 
surface capacity of gear systems depended on empirical formulas based on 
field experience. A major drawback in the application of the empirical 
methods is that they do not consider the effects of lubrication, which 
has been found to be of great importance in scuffing prediction. In 
recent years, the fully developed theory in elastohydrodynamic 
lubrication was successfully applied to spur and helical gears (9,52). 
Because of the geometric complexity of spiral bevel gears , an analysis 
covering all the possible factors to evaluate the risk of scuffing in 
lubricated spiral bevel gear drives has yet to be developed. 

In this chapter, an analysis of the temperature and film thickness 
for spiral bevel gear teeth is performed. The bulk temperature is calcu- 
lated by using the 3-D finite element method combined with a prescribed 
heat input. For the heavily loaded system, the limiting shear stress 
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depending on temperature and pressure is used to calculate the total 
heat generated over the whole contact ellipse. The average and maximum 
flash temperature in this contact ellipse is obtained by making a 
finite number of stripes. Each stripe is oriented along the direction 
of sliding velocity. It is assumed that the temperature profile along 
any stripe is the same as that of an infinitely long band heat source 
(in the direction perpendicular to sliding) whose width is equal to the 
stripe length and has the same heat flux profile along the stripe. The 
heat flux partition coefficient is calculated by assuming that the total 
heat is generated in the mid-plane of the lubricant film, and that there 
is no temperature jump in this plane. Hamrock-Dowson's equation is used 
to predict the film thickness between matching teeth. All quantities 
mentioned above, which include bulk temperature, flash temperature, 
traction coefficient and film thickness, are interdependent. They 
are solved as a system by an iterative method. 
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4.2 Flash Temperature 

During meshing, each tooth face will show a temperature increase 
(flash temperature) due to the frictional heat developed at the surface 
as the contact area moves along the tooth face. This temperature rise is 
restricted to the instantaneous contact area and will disappear very 
rapidly as soon as the area of tooth face is out of contact. Usually, 
this temperature is very high and it is an important factor in gear 
scuffing. 

■phe first successful prediction of flash temperature , done 
by Blok (13,14), was based on the heat conduction analysis of a 
semi- infinite body with a uniformly distributed moving heat source. 

J. C. Jaeger in 1942 (15) dealt with the problem of a moving source of 
heat with variable heat source and variable velocity. J. F. Archard 
(16) introduced a simple harmonic mean to obtain the interface temperature. 
A refined solution including a local heat partition coefficient between 
a pair of disks was derived by Cameron, et al. (17). More recently, 

Francis (18) made a further refinement in Blok's calculation by con- 
sidering a variable heat flux in the contact. 

The calculation of flash temperature is based on the assumption 
that the heat source passes over the surface of a semi- inf inite solid. 

For spiral bevel gears, the area of sliding contact is formed by elastic 
deformation of curved tooth surfaces where, according to the classic 
equation of Hertz (19), the area of contact is elliptical. Archard (16) 
has shown that, when the dimensionless parameter vR/a > vlO.O (where a 
is the thermal diffusivity, v is the sliding velocity and R is the 
nominal length), the heat flow in the direction perpendicular to sliding 
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may be neglected. The temperature distribution within a heat source 
of finite area can be determined by dividing the whole contact area 
into differential stripes parallel to the sliding direction. The 
temperature profile along any stripe is the same as that of an in- 
finitely long band heat source (in the direction perpendicular to slid- 
ing) of width equal to the stripe length and having the same heat flux 
profile along the stripe. 

The situation is illustrated in Fig. 4-1, which shows the 

contact heat flux distribution in an elliptical contact area with an 

angle a between direction of the velocity of heat flux (sliding velocity) 

and semi-minor axis of this area. Letting the equation (X/AMAX) 2 + 

2 

(Y/BMAX) = 1 describe this elliptic contact, the length J of a stripe con 
taining point (x,y) in the contact is 

YA = BMAX 2 tan 2 0 + AMAX 
YB = tan 9 BMAX 2 x^- - y_^ BMAX 2 tan 2 9 

YC = BMAX 2 tan 2 0 y 2 - 2 tan 0 BMAX 2 x ± y. 

+ BMAX 2 x 2 - AMAX 2 BMAX 2 

+ >- _ - YB ± VyB 2 - YA YC 

^ YA 

v' = tan 9 (y*’~ - y ± ) + x ± 


(4.1) 

(4.2) 

(4.3) 

(4.4) 


(4.5) 
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21 = ^(y* - y ± ) 2 + (x+ - x ± ) 2 (4.6) 

The center of this stripe is 

x. = 0.5 (xt + x ) (4.7) 

io i i 


y io = °' 5 (y i + V 


(4.8) 


The dimensionless parameter u is defined as 


u = ± -y(x- 


x io )2 + (y - y io )2 / 4 


(4.9) 


If (x,y) is in the forward position of this stripe (upstream of the 
sliding velocity), u is positive; otherwise, u is negative. The 
temperature at this point can be obtained by using the equation developed 
by Jaeger (15) for the temperature distribution along a stripe for a 

r* 

ifarger value of v£/a. 


T(x,y) 


la 

k 




a(H-u) 

irv 


(4.10) 


where k is the thermal conductivity and q is a constant heat flux 

(the reason for using a constant heat flux will be shown later) Cameron, 

et al. (17) showed that, for a uniform source, the temperature profile 

is given almost exactly for 10 < vJl/a < « by multiplying the asymptotic 

3/2 

(vi,/ a 'v °°) factor 1 + (2JL/i,-u)^*^[0.65 (a/v£) + 0.44 (a/vfc) ], and 


the temperature profile becomes 
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T(x,y) = 


2£ famuli 0 * 5 

k irv J 


r 2 ^ 0 - 5 o/o i 

1+ [^j [0.65 (a/v£) + 0.44(aM) 3/2 ]l 


(4.11) 


This factor will increase the maximum temperature and moves it toward the 
center of the stripes as vi/a decreases. Even if vl/a < 10, which occurs 
in the outer stripes of the contact ellipse, the equation can still be 
used because (1) the correction factor will hold approximately for 
vi/a < 10.0, and (2) the outer stripes are only a small fraction of the 
area of the heat source. Thus, there will be a minimal effect on the 
temperature in the whole contact area* 

Heat is generated by viscous shearing of the lubricant between the 
sliding surfaces. This heat is carried away either by the lubricant 
through convection or by the gear blanks through conduction. The 
relative importance between these two modes of heat transfer in EHD 
contact was examined by Trachman (20). He concluded that the heat 
carried away by heat convection can be neglected in comparison to the 
heat diffused to the gear body, except at very high speeds which exceed 
those of gears in current practice. Let the subscripts 1 and 2 denote 
the pinion and the gear, respectively. The heat generated in the contact 
area is divided betwen the pinion and the gear by the ratio a/(l - a), 
where a is the heat partition coefficient. Therefore, within the 
contact, the temperature distribution in the pinion and the gear are 


T l<*,y> = a 


KW] 

0.5 

L 

r 2& i ' 

TTV, 

< 

Z. — u. 

1 

V J 

\ 


1 1 


) 




0.5 


[0. 65 (a^v^, 


+ 0.44 a 1 /v 1 £ 1 ) 3/2 ] ( 


(4.12) 
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T 2 (x,y) 


= (1-ct) 


la 


a 2 ( ' l 2~ n l > 


trv. 


0.5 


1 + 


f U 2 l 0 - 5 


lV u 2 


[0.65 (a 2 /v 2 ^ 2 ) 


+ 0.44(a 2 /v 2 £ 2 ) 3/2 ]' 


(4.13) 


The lubricant local velocity profile and temperature profile 
for the very high sliding velocity case were estimated by Plint (21) 
and later confirmed by Trachman's analysis (20). For sliding contacts, 
their results showed a sharp S-shape velocity profile across the film 
with a large velocity gradient at the mid-plane. The temperature is at 
its nmirfinmn at the mid-plane and decreases almost linearly to both 
surfaces. Since shearing of the lubricant occurs mainly in the mid- 
plane, most of the heat is also generated in this layer. This 
results in a triangular temperature distribution across the film. 

According to Francis' analysis (18), if the bulk surface 
temperatures for the pinion and the gear, T^ and T^ 2 , are different 
(see Fig. 4-2), then the interfacial temperature can be expressed 


(a + Aa) 


Z n, q ' 

m 1 m 


+ T B1 ■ T ndd (x > y) 


= (1-a-Aa) 


t 2 + 


(h-Z )q 
in 


+ T 


B2 


(4.14) 


where Aa is the unknown function of (x,y) which expresses the local 
deviation from heat flux distribution aq and (1 - a)q, is the 
distance from pinion surface to the layers at which all the shear 


127 


is concentrated (Z^ - z >, k Q is the oil thermal conductivity, and 
T mid^ X ’ y ^ is the tem P era ture at the shear plane. Then 


(a + Aa) 


hq 


T + 

1 2k. 


+ T gl = (1 - a - Aa) 


hq 

X + — — 

2 2k 


+ T b2 ( 4 -15) 


A + Aa = 


T, + hq/2k + T - T 
t - o B2 B1 

T, + T, + hq/k 
1 ^ o 


(4.16) 


T mid <X ’ y) 


T 2 (x > y) + hq/2k 0 + T b2 - T B1 

T-^x.y) + T 2 (x,y) + hq/k Q 


( 

T x (x,y) 



+ T 


B1 


(4.17) 


The heat partition coefficient ct must now be determined to calcu- 
late the flash temperature of pinion and gear. The analytical expres- 
sion for a in terms of local conditions is difficult. The best approxima- 
tion of a will be the value that satisfies the condition for which the 
temperature fields of the two surfaces are identical over the contact 
ellipse. Thus, an analytical approximation can be written by equating 
the mean temperature of the two surfaces, as done by Jaeger (15) 


T 

l,ave 2k 


+ T B1 " - «> 


m Qh 

2,ave 2k 


+ T. 


B2 


(4.18) 


where 
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// T 1 (x,y)dxdy 

t = (4.19) 

l,ave it AMAX BMAX 

ff T,(x,y)dxdy 

T = 1 (4.20) 

i 2,ave tt AMAX BMAX 

Q is the total heat flux over the whole contact ellipse. Once 
the a is determined, the flash temperatures become 


T l f (x,y) = aT 1 (x,y) + < 4 * 21 > 

T 2 >f (x,y) = (1-a) X 2 (x,y) + T b2 (4.22) 
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4.3 Bulk Temperature 

The method mentioned in Section 4.2 is for calculating, with known 
heat flux, the flash temperature rise above the surface temperature which 
was assumed known beforehand. The existing criteria for scuffing of 
heavily loaded gears is that the total surface temperature rises to a 
critical value at which the oil film collapses, the protective surface 
layer breaks and metal to metal contact occurs. Uncertainty in the bulk 
surface temperature makes any scuffing prediction useless even if an 
accurate calculation of flash temperature can be calculated. 

Before the gear system starts to operate, all elements are at the 
ambient temperature. Once the gear system starts to operate, the 
temperature builds up as the gears are running due to frictional heat. 
After a sufficient period, the gear surface temperatures reach a steady 
state at which the heat flux flowing into the body is equal to that 
flowing out of the body. At each revolution, the tooth is subject to the 
same heating flux. Since the time period for each contact point in the 
contact zone is only a small fraction of the entire period of revolution, 
the local temperature jump (flash temperature) has completely decayed 
before it enters into the contact zone at the next revolution. Thus, it 
is justified to use an average heat input over the revolution to 
calculate the steady-state temperature rise of the body. 

The heat input is due to the heat generated at the instantaneous 
contact ellipse, and the amount depends on the load and the viscous 
shearing of the lubricant. The heat flux flowing out of the body is due 
to the heat convection to the surrounding air and lubricant. The 
relative importance of the heat transfer coefficient at different surface 
areas was discussed by Patir (22) and Townsend and Akin (23) in spur gear 
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systems. They also revealed the significant effect on temperature 
distribution by the lubrication method. In this study, the oil jet 
impingement depth is assumed to cover the whole area of contact side, 
which can be obtained by using the proper pressurized oil jet. Estimated 
values of the heat transfer coefficients are used on the other surface 
areas to calculate the bulk temperature. 

^ three-dimensional finite element program is used to calculate 
the temperature coefficient. The mesh of the system includes the gear 
shaft, gear body and the contact tooth with one adjacent tooth on both 
sides. The eight-node element is used. This model is the same as that 
used for the deflection coefficient except that the boundary conditions 
are different. In the temperature analysis, all the surfaces are subject 
to heat convection with different heat transfer coefficients except for 
the surface in the inner cylinder of the shaft, which is assumed to have 
no heat convection. The inner cylinder of the gear, created for the 
modeling, is very small in radius, therefore heat transfer can be 
neglected for these surfaces. Since the air in the inner cylinder of the 
pinion shaft is enclosed, there is no heat loss or gain in the steady 
state, and the heat convection at this surface can also be neglected. 

The heat transfer coefficient hj is assigned to the contacting tooth face 
^tiich is oil jet cooled. The topland and bottomland of the tooth and the 
other side of the tooth surface, which are not cooled by the oil jst, 
will have a heat transfer coefficient h t for air or air/oil mist. Since 
only three teeth are made in the model, there is a surface region A 
covering the surfaces where the teeth are taken off and the bottomland 
between those teeth. The heat transfer coefficient at this region A is 
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given a value of hj , which is the same as that of the surface coefficient 
cooled by the oil jet. The reason is that there is an oil jet cooled 
surface on each tooth and most of the heat will flow out of the tooth through 
this surface (h^ >> h fc ). All the other convective surfaces of the gear 
system are given a coefficient h g (Fig. 4.3). The theoretical estimated 
values of h g and tu can be found in (Refs. 24 and 25). However, the 
estimated values of h g , h and h fc based on the experimental results are 
used in this study (23). 

Thirty nodes are created on the contacting surface. The tempera- 
ture distribution on this surface is represented by T which is the 

JX 

temperature at the grid node J due to a unit heat flux at node I. By 

interpolation, the temperature at any contact point M due to a unit heat 

flux at the contact point N can be obtained in terms of T . The 

•J X 

details of this interpolation can be found in Appendix D. Once the 
contact path is located and the heat flux flowing into each body at each 
contact point has been calculated, the bulk temperature at the contact 
point M can be found as 

T M,B = z tC MN q N N=1...KMAX (4.23) 

where KMAX is the total number of contact points along the contact path. 
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4.4 Traction Coefficient 

In analyzing lubrication related failure mode in EHD contacts, the 
heat generated by traction in the contact area is the main concern. In a 
conventional analysis, the rheology of the lubricant is ignored and the 
local shear stress in the contact area is assumed to be proportional to 
the local pressure and acts in the direction opposite to the local slip. 

However, recent research into the shear behavior of EHD makes it possible 
to examine the role of the fluid rheology. 

The magnitude of the coefficient of friction in a highly loaded 
disk machine was examined by Misharin (26) at high rolling and sliding 
speeds. He also gave an approximated formula for the prediction of the 
friction coefficient based on his experimental data. But his formula did 
not give realistic results at a very high or very low value of the 
rolling or sliding speeds. Smith (27) presented a set of results that showed that 
the friction coefficient would increase to a maximum and then decrease 
as the sliding speed is increased. He pointed out that the decrease of the 
friction coefficient was due to thermal and non-Newtonian effects. 

The thermal effect upon the viscosity of lubricant was examined 
by Crook (28,29) and by Cheng (30). Crook used a disk machine to study 
rolling friction and sliding friction. He observed that the plot of 
traction against slip always showed a same kind of characteristic pattern 
as that which was found by Smith (27) regardless of the variation of load 
and speed. He also attempted to predict the traction analytically by 
simplified friction theory based on the following assumptions: The film 

thickness in the contact area is uniform; the pressure distribution in the 
contact area is Hertzian; the heat carried away by convection through the lubrican 
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is neglected; and the temperature rise on the surface is neglected. His 
method can predict a traction coefficient consistent with experimental 
data in high sliding region. However, he could not explain the behavior 
in low sliding region. Cheng used a thermal EHD theory which included 
energy equation and was free of all the assumptions made by Crook. 

Cheng’s results, compared with Crook's experimental data, still showed 
s big difference in friction coefficient in the low sliding region. Both 
suggested that the thermal effect alone cannot account for the 
sharp reduction of the effective viscosity in the low sliding region. 

^ ®f fort was made by Dyson (31) . He studied various 

experimental data and found that the curve of traction coefficient, when 
plotted against sliding speed, can be divided into three regions (Fig. 

. Region I is the linear region which shows the characteristic 
behavior of a Newtonian fluid. Region II is the nonlinear region which 
shows nonlinear variation of the traction coefficient with the shear rate, and 
region III is the thermal region which is dominated by the thermal effect 
on the shear stress of the lubricant. The fluid behavior in region I 
is shown to be explainable on the basis of linear viscoelasticity. In 
region II, the viscosity will decrease when the material undergoes a 
steady continuous shear with large strain. Dyson (31) and Gruber and 
Litvitz (32) suggested that the nonlinear behavior of the traction curve 
is due to this effect. At a high pressure, the relationship between 
stress and strain for the fluid is nonlinear; several equations were 
developed to describe this nonlinear behavior. Bell, Kannel and Allen 
(33), Hirst and Moore (34) used a "sinh" relation, and Trachman and Cheng 
(35) used a hyperbolic model. Observing that the traction force never 
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seems to exceed one tenth of the normal load. Smith (36), Johnson and 
Cameron (37), and Plint (38) proposed that the lubricant film would fail 
to retain a very high shear stress at extremely large strain rate, and 
shears like a plastic solid. 

More recent work on the rheological experiment was done by 
Johnson and Roberts (39). They observed that the liquid-solid transition 
was occurring in EHD contact and the transition point was dependent on 
temperature and pressure. The same observations were made by Johnson 
and Tevaarwerk (40) and Bair and Winer (41,42). Johnson and Tevaarwerk 
(43) then developed an elastic-perfect plastic solid model for the traction 
drive analysis. Based on their experimental results, most recently 
Bair and Winner (44) proposed the following model to predict the lubri- 
cant behavior over the whole range of visco-elastic-plastic change. 

4.4.1 Bair-Winner Model 

Bair and Winner fabricated an apparatus to measure the viscosity of 
the lubricant at high strain rates and a large shear stress similar to 
those in an EHD contact. They found that the lubricant under high pres- 
sure exhibits a classical visco-elastic behavior for small strain. For 
large strain and large stress, the lubricant behaves like a material 
having a limiting yield shear stress. For lubricant under a more 
moderate pressure but large strain, the result is that the shear stress is 
smaller and the lubricant is viscous. 

From their experimental results for 5P4E Santotrac-50 at large 
strain, the relationship of the dimensionless shear stress and the 
dimensionless shear rate can be described reasonably well with a single 
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natural log function. The dimensionless shear stress t is defined as 
the actual shear stress t(P,T) divided by the limiting shear stress 
t L (P ’ T >> an< * t * ie dimensionless shear rate y is the actual shear rate y 
multiplying by the viscosity p^CP.T) and divided by the limiting shear 

stress t^(P,T). t(P,T), t^(P,T) and y Q (P,T) are all pressure and 

temperature dependent. 

This natural log function can be expressed by 

y = - ta(l - x) (4 . 2 < 

which governs the nonlinear viscous flow. 

Based on the conventional Maxwell visco-elastic model, the total 
shear rate of the lubricant consists of the shear rate of elastic part 
(Y e ) and that of the viscous part (y ) 


Y = y + y 
'e 'v 


The elastic part can be nondimensionalized the s ame way 


; _ Y e » 0 (P,T) _ f(P,T) U q (P,T) ? 
Y “ i: l (P,T) " 6 o# (P,T)t (P,T) = T 


where G^CP.T) is the limiting elastic shear modulus at a given tempera- 
ture and pressure. With the above natural log relation for the viscous 
P ar t, a modified Maxwell model is obtained: 


Y = x - £n(l - t) 




T ' 

o L 
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In order to use this equation, three primary physical properties are 
required. These include the low shear stress viscosity y , the 
limiting yield shear stress t l> and the limiting elastic shear modulus 

G • These three parameters are all functions of the temperature and 

00 

pressure • 

In the application of the spiral bevel gears in heavy loaded 
cases, the pressure in most contact ellipses is very high (1 G Pa to 2 
G Pa or 150,000 psi to 300,000 psi) and the slide to roll ratio is about 
0.05 - 0.3. Under these severe conditions, the limiting shear stress 
for most lubricants should occur in most of the contact area. For the 
temperature calculation, it is safe to assume that the uniform limiting 
shear stress is applied at a given average pressure and temperature in 
the contact ellipse. 


T 


(P ,T ) 
ave ave 


(4.28) 


where t is given by experimental data. 

Li 

Since experimental data for t, in the temperatures encountered in spiral 

Li 

gears are not available, the conventional way to calculate this yield shear 
stress postulated by Dyson is used: 

T, - 6 > (4.29) 

L 00 

where G w is called limiting shear modulus, and is a function of 
temperature and pressure. Based on the experimental work of Hutton (11) 
and Switch (27), Dyson proposed a relationship between G^, temperature 
and pressure for high viscosity index mineral oils: 


137 


00 ’ T 2.52 + 0.0133 (T - 492) ' 1 * 45 x 1()4 (A. 30) 

where P is in psi, and T is in °R. 

A subroutine was prepared to evaluate the traction based on Bair 
and Winer’s model. Unfortunately, for the temperatures encountered in 

spiral bevel gears, the values of U(J , ^ and G. are not available from 
Bair and Winer ? s work. 
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4.5 Film Thickness 

In 1916, Martin provided the first paper about roller lubrication. 
Although he failed to predict a correct film thickness by using the 
assumptions of a rigid body and an incompressible isoviscous fluid, his 
approach was welcomed by later investigators. The influence on the film 
thickness by the pressure dependence of the viscosity in the contact area 
was examined by Gatcombe (45) by using an exponential relationship and by 
Hersey and Lowdenslager (46) by using parabolic viscosity-pressure 
relationships. More investigations in this field were made by Cameron 
(49) and by McEwen (48). These results showed that the effect of 
pressure dependence on viscosity could not alone account for the 
difference between prediced and observed film thickness, although it did 
provide some improvement. The effects of elastic deformation and 
pressure dependent viscosity were considered first by Ertel and Grubin in 
1945, who initiated the study of EHD lubrication. By assuming the 
deformation in the inlet region of contact zone was the same as that of 
the dry Hertzian contact, they calculated the pressure and film thickness 
in this region. Dowson and Higginson (19) developed a general formula 
for minimum film thickness for line contacts which included the effects 
of elastic distortion and pressure dependent viscosity. The dimension- 
less form of this equation is 


H 


min, l 


1.6 G 0,6 U 0,7 


W, 


0.13 


G = oE 

t 

U = n u/E R 

o 

w = w/e'r 

D 


(4.31) 


where 
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R is the effective radius of the roller pair, E' is the effective elastic 
modulus, n Q is the atmosphere viscosity of the lubricant, a is the 
pressure dependent viscosity parameter, n = n Q exp(ap), w is the load per 
unit length, and u is the speed. The thermal effect due to sliding 
was examined by Cheng and Stemlicht (50) by using numerical technique 
to solve the Reynold, elastic, energy and heat transfer equations. Their 
results indicated that the isothermal Dowson-Higginson equation could 
still predict a good approximation to the measured film thickness if n 

o 

was taken as the value at the bulk surface temperature of the contact 
bodies . 

The 'cylinder model' used in the lubrication analysis of gears 
is only good at the pitch point. Actually, the effective radius and rolling 
velocity are varied from point to point along the contact path. Wayne 
and Rodzimovsky (51) were among the first investigators who examined the 
actual contact ratio, actual involute profile, the combined effects of 
rolling and sliding, and film thickness in the whole course of the 
contact path. Unfortunately, their analysis did not include the elastic 
deformation in the contact zone. Based on the quasi-steady state assumption, 
Gu (52) used the EHD lubrication theory to calculate the film thickness 
of the spur gears by including the time dependent term 3(ph)/3t in the 
Reynolds equation. He found that this term did not have much effect on 
the minimum film thickness . 

In spur gears, the contact between mating teeth is a line 
contact extended from one end of the tooth to the other end of the tooth. 
Provided that the bulk surface temperature of the gear teeth is known 
beforehand, one cannot predict accurately the 


minimum film thickness by 
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using the Dows on -Higgins on equation. For conjugate spiral bevel gears, the 
contact is still a full line from one limit to the other limit. However, 
fully conjugated spiral bevel gears seldom exist in practice. The actual 
spiral bevel gears are crowned in both directions (lengthwise and profile- 
wise) to reduce sensitivity to the shift of contact patterns due to 
misalignment. In such cases, an elliptical contact area appears between 
the mating teeth under load, for which line contact solution is no longer 
valid . 

Hamrock and Dowson (53) evaluated numerically the analysis of an 
isothermal EHD lubrication point contact. They showed that the minimum 
film thickness could be related to the well-known line contact solution 
by simple expression involving the ellipticity parameter (54). In this 
study, Hamrock-Dowson' s point contact solution will be adapted in the 

elliptical contact of spiral bevel gears. 

The effective situation of contact between the spiral bevel gear and 

pinion can be seen in Fig. 4-5, in which there is a flat plane contact 

with a body which is described by the difference between neighboring 

surfaces of the gear and the pinion at the contact point. This curved body 

has effective radii R and R along the principal axis x and y, respec- 

x y 

tively . Under a load P, the surface near the flat plane will deform 

to an elliptical shape with semi-major axis AMAX and semi-minor axis 

BMAX. V and V are the velocities of the pinion and gear at the contact 
P S 

point. The ellipticity parameter was defined as the ratio of AMAX to 
BMAX. The minimum film thickness in the contact zone was related to 
Dowson-Higginson's line contact solution by the equation 
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H ^n- 


min min,£ 


where H 


min.il : Dimen sionless film thickness of Dowson-Higginson solution 


K : Ellipticity parameter 

H . ^min, £ 

min.Jl = — 

x 


U = V cos 0 
px p p 


U = V sin 6 

py P p 


U * V cos 0 
8 X g g 


D sy ' v g sln e g 


u - (U + U )/2.0 
x px gx" 


U = (U + U )/2 .0 

y py sy 


V =-^Ux 2 + Uy 2 

The dimensional form h ^ g of the minimum film thickness for a line contact can 


be written as 


h min,Jt " 1 * 6 «°‘ 6 (n o V) 0>7 e' 0,43 /w°* 13 


(4.33) 


where w P/2AMAX is the load per unit length and P is total load. The 
viscosity n Q is taken at the bulk surface temperature to account for the 
thermal effects on the film thickness. 
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It is important to note that n Q , P, AMAX, BMAX and V are varied 
along the path of contact. V depends on the gear kinematics, and AMAX 
and BMAX depend on the gear geometry and dynamic load P. n Q is strongly 
dependent upon the local static surface temperature which, in turn, is 
influenced by the local film thickness through frictional heating. Thus, 
the film thickness and the static (bulk) surface temperature are 
interdependent, and must be solved as a coupled system. The dynamic load 
is assumed to be uninfluenced by the film thickness. 
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4»6 Results of Lubrication Perf orman op 

The same set of gears used for the dynamic load calculations are 
used here to demonstrate calculation of the lubrication performance. Results 
were obtained for a range of operating conditions to determine the effect 
of speed, load, lubricant viscosity, and ambient temperature on the film 
thickness, bulk temperature, and flash temperature. 

The sliding velocity decreases from the beginning of the contact 
path where the gear tip contacts the pinion root, until the contact point 
is near the pitch point where the sliding velocity becomes zero. Then the 
sliding velocity Increases all the way to the end of the contact path 
where the pinion tip contacts the gear root. The current set of gears 
has the feature that the sliding velocity at the end of the contact path 
is larger than that at the beginning of the contact path; this fact 
creates a situation where more heat is generated at the end of the contact 
path. 

Typical distributions of the bulk temperature along the contact 
path are shown in Figs. 4-6 through Fig. 4-9 for various speeds. 

The bulk temperature of the pinion is always larger than that of the gear 
because the pinion speed is three times faster than the gear speed, and 
receives more heat in a unit time than the gear does. Although the 
temperature coefficients are higher near the gear tip, the maximum bulk 
temperature of both gears occurs at the end of the single tooth contact 
point where the maximum heat is generated. Distributions of the total 
flash temperature for the same cases are plotted on Fig. 4-10 through 
Fig. 4-13. The minimum flash temperature occurs at the pitch point where 
the sliding velocity is zero. For the high speed case, the variation 
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of dynamic load is less pronounced along the contact path. The rise of 
flash temperature on both sides of the pitch corresponds directly to the 
variation of sliding speed at the contact. The slight decrease at the 
end of the contact path is attributable to the decrease in dynamic load 
in this region. 

Fig. 4-14 shows the distributions of the film thickness for four 
different gear speeds. No excessive variations are seen along the 
contact path. A moderate peak is evident at the pitch point for the high 

speed cases, and this is associated with the slight drop of bulk tempera- 
ture at the pitch point. The steady rise of film thickness along the 
contact path is due to the increase in the entrainment velocity. The 
final uptrend of film thickness near the end of contact is again due to 
the decrease in the bulk temperature. 

Finally, the effects of an increase in ambient viscosity, ambient 

temperature, and load on the minimum film thickness maximum bulk 

and total flash temperatures is demonstrated in Fig. 4-15. As 
expected, an increase in viscosity would improve lubrication performance 
with a much greater film thickness and a slight drop in both bulk and 
flash temperature. An increase in ambient temperature would reduce 
the film thickness considerably, and an increase in load would also 
reduce the film thickness due to the increase in the bulk temperature. 
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"body 2 



4 —2 temperature distribution, across film 
thickness (h) with two different bulk 

temperatures 



high- slip 
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Fig* 4 -4 traction coefficient curve 
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Fig. 4-5. Simulated Contact Situation And Contact Area With 
Velocity of Pinion And Gear 
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Fig. 4-7 bulk temperature distribution for gear speed 150 
rad/sec 
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pig, 4-8 bulk temperature distribution for gear speed 300 
rad/sec 
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Fig. 4- 


9. Bulk Temperature Distribution For Gear Speed 
523 rad/sec 
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Fig. 4-10 flash temperature . distribution, for gear speed 50 
rad/ sec 
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Fig* 4 “11 flash temperature distribution for gear speed 150 
rad/ sec 



beginning of contact 


end of contact 


Fig. 4- 


12 flash temperature distribution for gear speed 300 
rad/sec 
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Fig. 4-14. Minimum Film Thickness is Plotted With Different Gear 
Speed 










Fig. 4-16 minimum film thickness, maximum flash temperature 
and maximun bulk temperautre for different 
ambient temperatures 
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CHAPTER V 
CONCLUDING REMARKS 

A computer solution for the dynamic load in a pair of spiral bevel 
gear sets was developed by solving the equations of motion for the pinion 
and gear shaft. An existing finite element code was used to calculate 
the combined stiffness of the contacting pinion and gear teeth as a 
function of contact position in the zone of action. In addition to 
the dynamic load analysis, a computer solution was developed to 
predict the bulk surface temperature, the flash temperature, and the film 
thickness along the contact path. An existing finite element heat code 
was used to calculate the temperature influence coefficients from 
which the bulk surface temperature is calculated. Both the lubricant 
film thickness and the sliding traction are calculated from the recent 
findings in EHL theories. 

Results were obtained for a set of experimental spiral bevel gears 
currently being tested at NASA Lewis Research Center. The results of 
dynamic load show that there exist numerous peaks in the variation of 
dynamic load against the gear shaft speed. These fluctuations correspond 
reasonably well with the critical frequencies of the system. The en- 
velope of the peaks suggests that the highest dynamic load occurs some- 
where near the critical frequency corresponding to the rotational mode 
oscillations of the two gears. 

Results of the film thickness show that its variation along the 
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contact path is not large, and is caused mainly by the increase in the 
entrainment velocity and the change in the bulk surface temperature. The 
total flash temperature variation is controlled by the sliding velocity, 
having its maximum near the end of the contact path where the transition 
from double to single teeth pair mesh occurs. Effects of operating 
variables on the minimum film thickness and maximum surface temperatures 
along the contact path can be obtained readily with this program. 

Results for the effect of change in the ambient viscosity show trends 
consistent with those anticipated from existing EHL theories. 
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TABLE II 

GEAR DATA AND LUBRICANT DATA 


Gear Data: 

Teeth 

Pitch Angle 
Shaft Angle 
Spiral Angle 
Diametral Pitch 


Gear Pinion 

36 12 

71°34' 18°26' 

90° 

35° 

5.14 


Standard Operating Conditions : 
Gear RPM 
Pinion RPM 

Load at Pitch Point N (lb) 
Ambient Temperature °C (°F) 


Geometric Dimensions: 
DGG = 0.1658 m 
ROG = 0.0312 m 
RIG = 0.04336 m 
RZG = 0.1964 m 
DGP = 0.2515 m 
ROP = 0.0325 m 
RIP = 0.09311 m 
RZP = 0.1987 m 


(see Fig. C-l) 
(6.527 in) 
(1.230 in) 
(1.707 in) 
(7.733 in) 
(9.901 in) 
(1.280 in) 
(3.6656 in) 
(7.824 in) 


5000 

15000 

11800 (2660) 
37.8 (100) 


Gear Material Data: 

Steel: ^ 3 

Density g/cm (lb/in ) 

Thermal Conductivity at 311°K (100°F) 
w/m°K (BTU/(pec)(in)(°F)): 

Young's Modulus GP fl (psi) 

Poisson Ratio 

Surface Convectivity . 

w/m 2 °K (BTU/ (sec) (in ) (°F)) : 

Oil Jet 
Oil/Air Mist 
Air 


7.81 (.282) 


46.7 (.000625) 
207. (30000000) 

0.3 


397. (0.000135) 

19.8 (0.00000675) 
3.97(0.00000135) 


Lubricant Data: 

Super-Ref ined, Napththenic, Mineral-Oil 
Dynamic Viscosity at 311°K (100°F) 
cp (lb* sec/in 2 ) 

Density at 311°K (100°F) g/cnr* (lb/in J ) 

Thermal Conductivity at 311°K (100°F) 

w/m°K (BTU/ (sec) (in) (°F)) 

Viscosity-Pressure • Temperature Relation 

v - U exp (ap + Bfcjr - ^r0) 

O L J-o 

Pressure-Viscosity Coefficient 
a m 2 /MN (in 2 /lb) 

Tempe r ature-Vis cosity Coefficient 
8 °K (°R) 


64.7 (0.00000938) 
0.61(0.022) 

0.125(0.00000168) 


0.023(0.00016) 
3890 (7000) 
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DC 


u 


d, e, c 


E 


m 


di 


NOMENCLATURE 
(For Chapter II) 

outer cone distance 

vector from machine center to the point on the surface 
addendum 

linear acceleration vector 

vector from wheel center to the point on the surface 
dedendum 

distance from cradle center to cam rotation center 
nominal wheel diameter 

machine coordinate 
side dresser offset 
blank offset 
running offset 
face width 

unit vector along the gear axis 

distance along tooth axis from crossing point 

inside dresser arm length 
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m. 

x 

N, n 


n 


n. 

x 


P 

Q 

R 

R 

Y 


Y 


d 


S 


(S) 

T 

T 

o 

*t 


u 


outside dresser arm length 
velocity ratio 

number of teeth for gear and pinion 
unit normal vector to the cutter blade 
index interval 

unit vector along the pinion axis 
cradle setting angle 
radius from tooth axis 

vector from crossing point to the point on the surface 

vector from the intersection of plane containing the 
cutter tips and wheel axis to the point on the surface 

side dresser radial 

vector from machine center to the intersection of plane 
containing the cutter tips and the wheel axis 

surface expression in matrix form 

cam setting 

standard cam setting 

unit vector along the direction of dressing diamond 
movement 

velocity of contact point 
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V 

v 

X, Y, X 
X 


*b 


di 


do 


X 


gy 


py 


6 

r, y 




r , y 

Y Y 


u 


linear velocity vector 

velocity of the contact point across the surface 

Cartesian coordinate 

machine center to back 

sliding base 

inside diamond setting 

outside diamond setting 

gear apex withdraw 

pinion apex withdraw 

side dresser axial 

cam guide angle 

eccentric angle 

pitch angle for gear and pinion 
face angle for gear and pinion 
root angle for gear and pinion 
cam rotation 
cam pitch radius 
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0 


0 


d 

X 



* 


wheel cutter rotation angle based on the line perpen- 
dicular to the S. 

dresser rotation angle 

work rotation 

radius 

shaft angle 

running shaft angle 

dresser block angle 

outside pressure angle 

cradle rotation 

spiral angle 

angular velocity vector 


Subscripts g and p refer to gear and pinion, respectively. 
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a 

F 

Fc 

Fr 


FMA 


±, j, k 
I 


J 


ks 


m 




RBG 


NOMENCLATURE 
(For Chapter III) 
acceleration vector 

external force vector 

tooth contact force 

bearing reaction force 

scalar quantity of contact force 

angular momentum about gravity center 

unit vector 

central mass moment of inertia about the principal axis 
of inertia 

equivalent mass 

combined tooth contact stiffness 
mass 

vector of external moment about gravity center 

position vector from gravity center to the tooth 
contact point 

position vector from gravity center to the bearing 
reaction point 

radius of effective base circle of gear 
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RBP 

T 

c 

x, y, z 

9 

n 


{ > 

[ 1 

[C] 

[DC] 

[DG] 

[DK] 

(G . 
ni 

{k} 

[«] 


(6 



radius of effective base circle of pinion 
static force 

translational displacement 
rotational displacement 

relative rotational velocity of moving coordinate to 
fixed coordinate 

nominal shaft velocity 

column matrix 

matrix form 

damping matrix 

transformation matrix 

transformation matrix 

transformation matrix 

normal displacement at point i due to rigid body motion 
stiffness matrix 
mass matrix 

normal displacement at point i due to tooth contact 
force 



177 


transformation matrix 


transformation matrix 


Subscripts p , g refer to pinion and gear. 
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a 


AMAX 

UMAX 

E' 

G 

G 

00 

h 

h j 

h 

s 


H . 
mm 

^min , l 
k 

k 

o 

K 

i 


NOMENCLATURE 
(For Chapter IV) 

thermal diffusivity 

length of semi-major axis of contact ellipse 

length of semi-minor axis of contact ellipse 

effective elastic modulus 

material parameter 

limiting elastic shear modulus 

film thickness 

surface heat transfer coefficient for oil jet cooling 

surface heat transfer coefficient for free convective 
cooling 

surface heat transfer coefficient for air/oil mist 
cooling 

dimensionless film thickness of point contact solution 

dimensionless film thickness of line contact solution 

solid material thermal conductivity 

oil thermal conductivity 

ellipticity parameter 

length of contact stripe 
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P 

q 

Q 


R 


ave 


T 

mid 


u, V 


U 


W„ 


m 


a 


pressure 
heat flux 

heat flux at tooth contact point 
nominal length, effective radius 
temperature 

average temperature over whole contact ellipse 

bulk temperature 

flash temperature 

temperature at shear plane 

speed 

speed parameter 
sliding velocity 
load per unit length 
load parameter 

distance from solid surface to shear plane 

heat partition coefficient, pressure dependent viscosity 
parameter 

shear rate 
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Y 


Y 


Y 

n 

y 

T 

T 



dimensionless shear rate (y y o /T L ) 
shear rate of elastic part 
shear rate of viscous part 
atmosphere viscosity of lubricant 
low shear stress viscosity 
shear stress 

dimensionless shear stress (t/t l ) 
limiting shear stress 
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appendix a 

SAMPLE DATA FOR MACHINE SETTINGS 

Gear Pinion 




Concave 

Convex 

Number of Teeth 

36 


12 

Face Width 


1.00 


Shaft Angle 


1.5708 


Outer Cone Distance 


3.691 


Addendum 

0.093 


0.231 

Dedendum 

0.267 


0.129 

Pitch Angle 

1.249 


0.3217 

Face Angle 

1.276 


0.3895 

Root Angle 

1.1813 


0.2944 

Spiral Angle 


0.61086 


Nominal Wheel Diameter 


6.0 


Machine Center to Back 

0.0 

-0.067 

0.060 

Sliding Base 

-0.016 

-0.009 

-0.046 

Blank Offset 

0.000 

0.041 

-0.045 

Cam Setting 

6.916 

7.089 

7.279 

Eccentric Angle 

0.3598 

0.3616 

0.3604 

Cradle Angle 

1.2116 

5.3354 

5.4990 

Cam Guide Angle 

0.0 

0.0 

0.0 

Standard Cam Setting 

7.00 

7.00 

7.00 

Cam Pitch Radius 

6.2273 

5.9567 

5.9567 

Index Interval 

11 

11 

11 
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APPENDIX A (continued) 


Dresser Block Angle 
Outside Pressure Angle 
Side Dresser Radial 
Outside Dresser Arm Length 
Outside Diamond Setting 
Inside Dresser Arm Length 
Inside Diamond Setting 
Side Dresser Axial 
Side Dresser Offset 


Gear Pinion 



Concave 

Convex 

0.7854 

0.7854 


0.3490 

0.3490 

0.3482 

3.113 

3.028 

2.959 

1.0 

1.0 

1.0 

0.0 

0.0 

0.171 

1.0 

1.0 

1.0 

0.223 

0.170 

0.0 

0.0 

0.023 

0.029 

1.00 

0.800 

1.200 
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APPENDIX B 

From the geometry of Fig. 2-10, the motion formulas for 
3rd cam rotation are obtained as follows! 

, n u 

*0 dc + n e o 


An 

sin 1 K - sin ip + sin a + — £ sin(e - ip ) 
^ 0 DC O O 


* - * 2 ~ 2 

e 2 = e o + *2 + *0 


e = ^2 ~ a 


n i 

A = — e 
n 


From (B-l) 


dij/ n 

o u 

de DC + n 
o u 


Differentiating (B-2) with respect to s q 


d’i'2 An di(> 

cos ^9 TT~ = cos 'P T"^ + cos (e -ip ) (1 - -r-2-) 
4 de o de DC o o de } 

o o O 


ife ' 3T s “ *2 


u 


DC + ti 

u 


cos \p + 


An 


„ nr ,“T cos(e - ip ) 

o DC + n o o 


u 


sec ip 2 
DC + n 

u 


stand- 

(B-l) 

(B-2) 

(B-3) 

(B-4) 

(B-5) 

(B-6) 


= N 


(B-7) 
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where N = n u cos + An u cos(e Q - i|> Q ) 


From (B-4) and (B-5) 


^ = . i + 1!* . !!° 

de de de de 

o o o o 


sec ip r 


= 1 + N 


DC + n DC + n 
u u 


= D 


sec 

i 

DC + n 


u 


(B-8) 


dX 


where D = DC cos i|> + N and also ^ 


n. 

l 

n 


de 

de. 


, de , , de 

dip _ dip o _ djp o n 

dX ~ de dX de de n. 

o o i 


_n_ N 
n i D 


(B-9) 


The term a = 
c 


f ^ 

a 2 ii) 

,3X 2 , 


used In determining the motion parameters is 


derived as follows: 


d 2 t{i 

dX 2 


n 1 L dN dD) 

l D dx ' « dxj 



2 , 

n 

1 

n 

2 

X 

V > 

D l 


d 35 . „ 

de de 
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)2 (DC+n ) 


cos 0 


de de 

o o 


and 


dN 

de 


dip 

/■ n 

r \ 

dip 

= ~ sin ^ ~a " An sin 

u 0 de u 

n 

e - ip 
0 0 

l s? 

U 

1 i 

1 0) 


AH 

DC + n 

u 


where 


AH = 


2 

Tl 

U 


sin ip + An DC sin(e - ip ) 
° o o 


- - DC sin *, “ 


2 de ~ DC + n 
o u 


5 c + n [DC N tan ^2 + m) 

U 


dh , n 2 cos *9 

— 9 = 5 — ^ (N(DC N tan ip 0 

a \ 1 n .- 2 


d\~ “i D 


(B-10) 


(B-ll) 


(B-12) 


(B-12) 


+ AM) - D AM} 


(B-13) 
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APPENDIX C 


The details of the mass, viscous and stiffness matrices of the 
equation of motion are shown here. 


MAXC 

AXYZCG (I,K) 

AXYZCP (I,K) 
AI (J,K) 

XGG, YGG, ZGG 

AXG, AYG, AZG 
SGP, YGP, ZGP 

AXP, AYP, AZP 
SYZCG (I,K) 

XYZCP (I,K) 


Number of contact points at each time 

Normal vector of gear surface at contact point 
(inward to gear body) 

Normal vector of pinion surface at contact 
point (outward to pinion body) 

Stiffness of contact point I due to a unit 
load applied at contact point K 

Translational displacements of gear gravity 
center 

Rotational displacements of gear body 

Translational displacements of pinion gravity 
center 

Rotational displacements of pinion body 

Position vector from gear gravity center to 
contact point K 

Position vector from pinion gravity center 
to contact point K 


where 1=1, 2, 3 indicates the component in three-dimensional axis. 

J and K equal one to MAXC. Array ABG (J,I), J = 1 to 12, I - 1 to MAXC 

was introduced for convenience in developing equations: 

ABG (1,1) = Z - AI(I,K)*AXYZCG(1,K) 

K 
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ABG (2,1) 

= Z - AI (I , K) *AXYZCG (2 , K) 
K 

ABG (3,1) 

= Z - AI ( I , K) *AXYZCG (3 , K) 
K 

ABG (4,1) 

= Z - AI(I,K) [XYZCG(2,K)*AXYZCG(3,K) 
K 


- XYZCG(3 ,K)*AXYZCG(2 ,K) ] 

ABG (5,1) 

= Z - AI ( I , K) [AYZCG (3 , K) *AXYZCG (1 , K) 
K 


- XYZCG (1 , K) *AXYZCG (3 , K) ] 

ABG(6,K) 

■ Z - AI(I,K)[XYZCG(1,K)*AXYZCG(2,K) 
K 


- XYZCG (2 , K)*AXYZCG (1 , K) ] 

ABG (7,1) 

■ Z AI (I , K) *AXYZCP (1 , K) 
K 

ABG (8,1) 

= Z AI ( I , K) *AXYZCP (2 , K) 
K 

ABG (9,1) 

= Z AI ( I , K) *AXYZCP (3 , K) 
K 

ABG (10,1) 

= Z AI(I,K)[XYZCP(2,K)*AXYZCP(3,K) 
K 


- XYZCP (3 , K) *AXYZCP (2 , K) ] 

ABG (11,1) 

= Z AI(I,K)[XYZCP(3,K)*AXYZCP(1,K) 
K 


- XYZCP (1 ,K)*AXYZCP (3 ,K) ] 

ABG (12,1) 

■ Z AI(I,K)[XYZCP(1,K)*AXYZCP(2,K) 
K 


- XYZCP (2 ,K)*AXYZCP (1 ,K) ] 

The summation is from K = 1 to K = KMAXC, and I = 1 to KMAXC. 

The force situation is shown in Fig. 3-1 and Fig. C-l. 
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For the gear, the components of the contact force are 


AXYZCG (1,1) {ABG(J,I)) T 



+AXYZCG 


T 

(2,1) (ABG(J , I) ) 


D I 

P' 


+AXYZCG 


(3,1) (ABG(J,I)} T 



I D 

v p 


J = 1 to 12 


where { 


denotes column matrix. 



= (XGG, YGG, ZGG, AXG, AYG, AZG, 




XGP, YGP, ZGP, AXP, AYP, AZP) and I from one to MAXC. The damping force 
at contact point I is 

the force terms due to bearing reaction are 


-2KXP*XGP + KXP(RIP-R2P)*AZP i 
-KYP*YGP j p 

-2KZP*ZGP + KZP (R2P-R1P ) *AXP k 

p 


and the torque due to these bearing forces are 
- (R1P-YGP o ) KZP*ZGP- (R 1P-YGP 0 ) 2 KZP*AXP 
+ (R2P+YGP ) KZP- (R2P+YGP 0 ) KZP + AXP i p 
+(RlP-YGPo )KXP*XGP- (RIP- YGP o ) 2 KXP*AZP 
= (R2P+YGPo ) KXP*XGP- (R2P+YGP <» ) 2 KXP*AZP k^ 
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where YGP 0 is the rigid body displacement of the pinion in the pinion 
fixed coordinate. 

The damping force due to the bearing is 

C' *XGP + C' *AZP T 
xp azp p 

+C' *YGP T 

yp j p 

• a 

+C' *ZGP + C’ *AXP k 
zp axp p 

and the damping torque is 

C" *ZGP + C" *AXP 7 
zp axp p 

+C" *XGP + C" *AZP k 
xp azp p 

Now the coefficients of damping will be determined. Let 1=1 

to MAXC and introduce array CBG(j,I), J = 1 to 12 for convenience, and 
define 

CBG (1,1) = -AXYZCG(1,I) 

CBG (2,1) = -AXYZCG(2 , I) 

CBG (3,1) = -AXYZCG(3 , I) 

CBG (4,1) = -[XYZCG(2, I)*AXYZCG(3 , I)-XYZCG(3 , I)*AXYZCG(2, I) ] 
CBG (5,1) = -[XYZCG(3,I)*AXYZCG(1,I)-XYZCG(1, I)*AXYZCG(3 , I) ] 
CBG (6,1) = -[XYZCG (1,I)*AXYZCG(2,I)-XYZCG(2, I)*AXYZCG(1,I) ] 
CBG (7,1) = AXYZCP(1, I) 

CBG (8,1) = AXYZCP(2,I) 

CBG (9,1) = AXYZCP(3,I) 

CBG (10, I) = XYZCP (2,1) *AXYZCP (3,1) -XYZCP (3,1) *AXYZCP (2,1) 
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CBG (11,1) = XYZCP (3 , I) *AXYZCP (1 , I) -XYZCP (1,1) *AXYZCP (3,1) 
CBG (12,1) = XYZCP (1 , I)*AXYZCP (2 , I) -XYZCP (2 , I) *AXYZCP (1,1) 


For gear, J = 1 to 12 

C = AXYZCG(1,I) (CGB(J,I)} T C 

cxg c 

C = AXYZCG(2,I) (CBG(J,I)} T C 

cyg 

C = AXYZCG(3 ,1) (CBG(j,I)} T C 

czg C 


C’ = [XYZCG(2,I)*AXYZCG(3,I)-XYZCG(3 ,1)*AXYZCG(2,I)] 
cxg 

{CBG(J,I)} T C c 


c’ = [XYZCG (3,1) *AXYZCG (1,1) -XYZCG (1,1) *AXYZCG (3,1)] 
cyg 

{CBG(J,I)} T C c 


C' = [XYZCG(1,I)*AXYZCG(2,I)-XYZCG(2,I)*AXYZCG(1,I) ] 
czg 

(cbg(j,i)} t c c 


C' = -COG 
xg 


C' = -C1G - C2G 

yg 


C' = C1G* (RIG-XGGo )-C2G* (R2G-XGG 0 ) 
azg 


c » = -C1G-C2G 
zg 


C' » -C1G* (R1G4XGG a ) +C2G* (R2G-XGG 0 ) 

ayg 


C" = -(R1G+XGG 0 )C1G+(R2G-XGG 0 )C2G 


C" = -(R1G+XGGJ 2 C1G- (R2G-XGG 0 ) 2 C2G 
azg 
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C" = (RlG+XGG o )ClG-(R2G-XGG 0 )C2G 

J © 

C ayg = -(^G+XGGo) 2 C1G-(R2G-XGG 0 ) 2 C2G 

For pinion, J = 1 to 12 

C - -AXYZCP(1,I){CBG(J,I)} T c 

exp C 

C = -AXYZCP(2,I){CBG(J, I) } T C 

r C 

C * -AXYZCP(3,I){CBG(J,I) } T C 
Cd P c 



The torque due to this contact force is 

[XYZCG(2,I)*AXYZCG(3, I)-XYZCG(3, I)*AXYZCG(2, I)] {ABG(J,I) T 


i 



+ [XYZCG(2,I)*AXYZCG(3 , I)-XYZXG(3 , I)*AXYZCG(2, I) ] {ABG(J,I)} T 



+ [XYZCG(1,I)*AXYZCG(2, I)-XYZCG(2, I)*AXYZCG(1, I) ] {ABG(J,I)} T 



and the damping torque is 


c-3 
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C' 

cxg 


• 


• > — 


D 

i + C * 

D t + C' 

D 

gl 

, g cyg , 

gig czg 

* 

D 


l D 1 

D 

Pi 


l p/ 

l P/ 


k 

g 


The force terms due to the bearing reaction are 

-KXG + XGG T 
g 

-2KYG*YGG+KYG* (R±G+XGG ) -KYG* (R2G-XGG) J 

6 

- 2KZG*ZGG-KZG* (R1G4XGG) *AYG+KZG* (R2G-XGG ) * AYG k 

g 

The torque terms due to these bearing reactions are 
-R1G+XGG 0 ) *KZG ZGG- (R1G+XGG 0 ) 2 *KZG*AYG 

+R2G-XGG 0 ) *KZG*ZGG- (R2G-XGG 0 ) 2 *KZG*AYG J 

o 

+(RlG+XGG 0 )*KYG*YGG-(R2G-XGGo) 2 *KYG AZG k 

o 

where XGG 0 is the rigid body displacement of gear in the gear fixed 
coordinate. 

The damping force is 


c 

xg 

c f 

yg 

c* 

Zg 


* XGG i 

g 


* ZGG + C ' 

azg 


* ZGG + C' 

ayg 


• 

* AZG j 

J g 

* AYG k 

g 


and the damping torque is 


c" * 

ZGG 

+ 

C" 

Zg 



ayg 

C" * 

YGG 

+ 

c" 

yg 



ayg 


* AYG i 

J g 

* AZG k 

g 
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For pinion case, the components of contact force are 


-AXYZCP(1,I; (ABG(J,I)} T ( 


D 

l P 


-AXYZCP(2, I) (ABG(J,I)} 


Pi 


-AXYZCP(3,I) (ABG(J,I)r 


v P/ 


the damping force at contact point I is 
C 


I = 1 to MAXC and 
J - 1 to 12 


fn 1 

i +C . • , 

j +c 

• 

8 

P cyp |D g | 

J P CZ P ! 

f D ) 

! g 

1 • 

Id J 

i D J 

| 

D 

p 

' p' 


1 P' 


kp 


the torque due to this contact force is 

- [XYZCP (2,1) *AXYZCP (3 , I)-XYZCP (3 , I) *AXYZCP (2,1) {ABG(J,I)} T 

I 

I 

g 

) 

P 

- [XYZCP (3 , I) *AXYZCP (1 , I)-XYZCP (1 , I) *AXYZCP (3,1)] { ABG ( J , I) J 1 
D 

8 
D 

l Pi 

- [XYZCP (1 , I)*AXYZCP (2 , I)-XYZCP (2 , I)*AXYZCP (1 , I) ] { ABG ( J , I) J 1 


'd ' 

g 


v P/ 
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J = 1 to 12 and the damping torque is 


+ C' 
cyp 


C' = -[XYZCP(2,I)*AXYZCP(3 ,I)-XYZCP(3,I)*AXYZCP(2,I) ] 
exp 

{cbg(j,i)} t c c 


C' = - [XYZCP (3 , I)*AXYZCP (1 , I)-XYZCP (1 , I) *AXYZCP (3,1)] 
cyp 

{ cbg(j,i)} T C c 


C' = -[XYZCP(1,I)*AXYZCP(2,I)-XYZCP(2,I)*AXYZCP(1,I)] 
czp 

{cbg(j,i)} T c c 


c' = -C1P - C2P 
xp 


c' = C1P* (R1P-YGP „ ) - C2P* (R2P+YGP 0 ) 
azp 


C' = -COP 

yp 


C' = -C1P - C2P 
zp 


C' = -C1P* (RIP-YGPo ) + C2P* (R2P+YGP 0 ) 
axp 


C" - - (R1P-YGP o ) C1P + (R2P+YGP 0 ) C2P 
zp 

C" * - (R1P-YGP 8 ) 2 C1P-(R2P+YGP 0 ) 2 C2P 
axp 


C" = (R1P-YGP o ) C1P- (R2P+YGP ) C2P 
xp u 

C" = - (R1P-YGP 0 ) 2 C1P- (R2P+YGP 0 ) 2 C2P 
azp 



j +C' 

J p czp 


p' 
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Pig. C-l 


Bearing Force Configuration 
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APPENDIX D 

INTERPOLATION OF THE COEFFICIENTS OF DEFLECTION 
AND BULK TEMPERATURE 

(A) Deflection Coefficient 

There are thirty nodal points in the contacting tooth surface. 

The deflection coefficient D is the deflection (normal direction) at 

point i due to a unit load applied at point j. The coefficient D^ was 

obtained from the finite element method. Fig. D-l shows these grid points 

in the axial plane and also shows the contact path. It is desired to 

r 

determine the deflection D at contact point I in surface element a 
with four corner points jl, j2, j3 and j4, due to a unit load at this 
point I. The concept of the shape function used in the finite element method 
was borrowed here to solve the problem* A mapping diagram is shown in 
Fig. D-2, correspondent to Fig. D-l. All the surface elements in 
this mapping diagram are of the same size as the square shape. Intro- 
ducing two new local coordinate variables r and s, one can obtain the 
deflection at I by 

Dj - 0.25*(l-r)*(l-s)*D J1 J;L + 0.25*(l+r)*(l-s)*D J2>J2 

+ 0.25* (1+r ) * (1+s ) *D j 3 J3 + 0.25*(l-r)*(l+s)*D j4 ^ j4 
r, s are defined as 

Xj = 0.25*(l-r)*(l-s)*X J1 + 0.25* (1+r )*(l-s)*X J2 


+ 0.25* (1+r ) * (1+s) *X J3 + 0.25*(l-r)*(l+s)*X j4 
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Y j " 0 • 25* (1-r ) * (1-s ) *Yj 1 + 0.25*(l+r)*(l-s)*Y J2 
+ 0.25* (1+r ) * (1+s ) *Y J3 + 0.25*(l-r)*(l+ s )*Y 

l ^ 

I 

which X and Y are the coordinates in the axial plane. 

(B) Bulk Temperature 

The temperature rise (tJ) at contact point I due to a unit heat 
flux at the same point can be calculated the same way as that for the 

deflection coefficient. If T„ is the temperature rise at nodal point i 
due to a unit heat flux at nodal point j , then 

Tj = 0.25*(1- )*(l-s)*T JljJ1 + 0.25* (1+r )*(l-s)*T J2>J2 

+ 0.25* (1+r)* (1+s )T + 0.25*(l-r)*(l-s)*T 

J4,J4 

Now the temperature rise (T® ) at contact point J due to a unit 

heat flux at contact point I (I^J) has to be determined. If point J is 

I 

far away from point I, then the temperature rise at point J will not be 
influenced by the local temperature jump at point I. Fig. D-3 shows 
this situation where contact I is in surface element g with nodal points 
II, 12, 13, 14 and contact point J is in surface element a with nodal 
point Jl, J2, J3 and J4. The temperature rise at nodal points II, 12, 

13, and 14 due to heat flux at contact point J is 

T II,J * °- 25 * <1 - r J ) * (1 - s j ) * T n,ji + 0.25*(l-tr J )*(l-, J )n IIfJ j 

' + °- 25 * (1+r J ) *< 1+s J>* T II,j3 + «.2J*(l-r J )*(14, J )*T IIfJ4 

I 

where II = II, 12, 13 and 14, respectively, and r , s T are local 

l J J 

i 

I 

I 

! 

i 
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coordinates in surface element a. Then the temperature rise at contact 
point I will be 

T I J = J + 0.25*(l+r I )*(l-s I )*T I2 ^ 

+ 0.25*(l+r I )*(l+s I )*T I3 ^ J + 0.25*(l-r 1 )*(l+s I )*T l4>J 

where y and S^. are local coordinates in surface element B. 

If contact point I is close to contact point J, the local 

temperature effect should be counted. One can make a square with the 

same size as any other squares so that the contact point X will be 

within this square and contact point J will become one of the four corner 

points of this square (see Fig. D-4). Then the temperature at each 

corner point due to heat flux at point J can be determined as above. 

Let r and s be the local coordinates in the new surface element, then 
I I 

T C T = 0.25*(l-r )*(l-s )*T + 0.25*(l+r )*(l-s )*T 

I, J J- 1J J-- 1 - 

+ 0.25*(l+r l )*(l+s I )*T 2 + 0.25*(l-r ;[ )*(l+s I )*T 3 

the order of Tj, T^, T 2 and T 3 will be changed, depending on how this 
new element is formed. In order to cover all the possible cases where 
this new element may locate, a new element layer was created to surround 
the original mapping diagram (as shown in Fig. D-4). The temperature at 
these new added grid points due to heat flux at any possible grid point 
was given artificially. 
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Fig. D-4. Mapping Diagram With A New Element Layer 
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